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ABSTRACT 


The  specific  DC  flow  resistance,  characteristic  impedance  and  complex  wave 
number  of  air-filled  fiberglass  have  been  extracted  from  measurements  of  the 
peak  and  half-power  frequencies  of  the  two  lowest- frequency  normal  modes  of  a 
rectangular  plexiglass  cavity  partially  filled  with  fiberglass.  Using  Biot 
theory  [Ref's  8,  9  and  10]  to  describe  the  sound  field  within  the  fiberglass, 
a  dispersion  relation  for  the  allowed  modes  in  the  cavity  was  derived,  from 
which  the  above  properties  were  extracted.  The  values  determined  by  this  new 
method  compared  favorably  (within  25%)  with  the  results  of  measurements  using 
classical  methods.  The  new  method  differs  from  the  classical  resonance  tube 
method  [Ref.  4]  by  1)  the  use  of  a  wave  theory  to  descibe  the  sound  field 
within  the  fiberglass  and  2)  the  use  of  a  mode  with  particle  velocity  parallel 
to  the  air-fiberglass  interface.  The  differences  make  it  possible  to  extract 
the  acoustic  properties  of  bulk  fiberglass  using  only  a  thin  sample. 
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I.   INTRODUCTION 

A.  BACKGROUND 

Elastic  wave  propagation  in  a  fluid-saturated  porous  solid  is  a  very 
important  and  active  subject  of  research  with  a  wide  spectrum  of  applications 
including  oil  prospecting,  nuclear  test  ban  treaty  verification  and  sound 
transmission  and  reflection  at  the  ocean  floor. 

When  the  fluid  in  the  pores  is  highly  compressible  (air  for  example)  the 
acoustic  properties  are  completely  determined  by  the  characteristic  impedance, 
zchar'  or'  equivalently  by  the  complex  wave  number,  k  .  From  these  guantities 
can  be  extracted  parameters  which  depend  only  on  the  microscopic  geometry  of 
the  porous  solid,  the  values  of  which  can  then  be  used  when  the  fluid  in  the 
pores  is  not  highly  compressible  (water  for  example).  It  is  therefore 
desirable  to  have  accurate  and  reliable  methods  by  which  these  parameters  can 
be  measured. 

B.  HISTORY  OF  CHARACTERISTIC  IMPEDANCE  AND  COMPLEX  WAVE  NUMBER  MEASURING 
TECHNIQUES 

Rayleigh  (1929)  [Ref.  1,  p.  328]  has  investigated  the  performance  of  an 

idealized  absorber  of  sound  consisting  of  small  bundles  of  tubes.  He  extended 

the  theories  of  Helmholtz  and  Kirchoff  that  relate  to  the  transmission  of 

sound  along  tubes  and  was  able  to  deduce  a  value  for  the  absorption 

coefficient  at  the  surface  of  such  a  material  in  terms  of  the  propagation  of 

the  surface  area  occupied  buy  the  open  pores,  the  radius  of  the  tubes,  and  the 

physical  constants  of  air.   This  analysis  has  provided  a  basis  on  which  to 

found  more  general  theories  and  description  of  sound  propagation  through 

porous  materials. 
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Scott  (1945)  [Ref.  2]  proposed  a  method  of  determininq  the  complex  wave 
number,  k  which  describes  the  spatial  variation  of  the  acoustic  disturbance  in 
a  porous  medium.  This  nethod  consists  of  measurinq  the  attenuation  constant 
and  the  wavelenqth  constant  of  a  plane  proqressive  wave  by  orobinq  the  porous 
medium  with  a  probe  tube  microphone.  Scott  then  determined  the  chatacteristic 
impedance  of  the  material  with  an  apparatus  devised  by  H.  0.  Taylor  (1913) 
[Ref.  3]  to  measure  the  sound  dampinq  effects  of  acoustical  materials.  This 
apparatus  is  similiar  to  the  present  day  impedance  tube. 

Beranek  (1949)  [Ref.  4]  describes  several  methods  of  obtaininq  a 
material's  acoustic  impedance  bv  measurinq  the  properties  of  the  sound  field 
in  the  vicinity  of  an  absorbinq  material  exposed  to  a  normally  incident  sound 
wave.  Two  of  the  most  tractable  methods  are  based  on  transmission  line 
theory:  the  impedance  tube  method  and  the  resonance  method  (normal 
incidence),  with  the  latter  method  qenerally  beinq  viable  only  for  determininq 
surface  impedance.  Beranek  also  developed  a  method  to  calculate  the  complex 
prooaqation  constant  from  measurements  of  a  material's  characteristic 
properties,  i.e.  porosity,  pore  size,  etc.  Once  these  properties  are 
determined,  the  characteristic  impedance  can  be  calculated. 

More  recently,  Yaniv  (1973)  [Ref.  5]  developed  a  new  method  of 
calculatinq  both  the  complex  propaqation  constant  y  =  jk  and  characteristic 
impedance  of  a  material  usinq  transmission  line  theory  applied  to  an  impedance 
tube  device.  His  method  consists  of  measurinq  standinq  wave  properties  in  the 
proximity  of  an  absorbinq  material  alternately  backed  with  a  hiqh  impedance 
and  a  low  impedance  termination. 

Finally,  in  1980,  Chunq  and  Blaser  [Ref.  6]  published  their  article 
discussinq  a  Dual-flicrophone  Transfer  Function  technique  for  measurinq  in-duct 
acoustic  properties,     namely   the  acoustic  impendance   and  reflection 
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coefficient  from  which  the  complex  wave  number  can  be  determined.  This 
technique  was  compared  to  the  conventional  Standing  Wave  Tube  technique  and 
found  to  qive  similar  results  [Ref .  7] . 

C.   PURPOSE  OF  STUDY 

The  objective  of  this  research  is  to  investiqate  the  feasibility  of 
obtaininq  the  characteristic  impedance  and  the  complex  wave  number  of  a  fluid- 
filled  porous  solid  by  measurinq  the  resonance  frequencies  and  half-power 
widths  of  a  fluid-filled  rectanqular  cavity  which  is  partially  occupied  by  a 
porous  solid  of  varying  thickness.  In  this  investiqation,  the  fluid  is  air, 
the  porous  solid  is  fiberqlass  and  the  cavity  is  made  of  plexiqlass. 

This  investiqation  differs  from  previous  invest iqat ions  in  that  1)  Biot 
Theory  [Ref.  8,  9,  10,  and  11]  is  applied  to  solve  for  the  allowed  modes  of 
the  partially-filled  cavity,  and  2)  for  the  primary  mode  of  interest  the 
particle  velocity  is  directed  alonq  the  lenqth  of  the  fluid-porous  solid 
boundary  as  opposed  to  normal  to  this  boundary. 

Table  1  summarizes  the  advantaqes  in  the  course  of  this  research.  The 
advantaqe  of  this  technique  over  the  currently  existinq  methods  is  that  it 
does  not  require  probinq  the  sound  field  either  within  or  in  front  of  the 
porous  solid  as  do  methods  1,2,  3,  and  5,  nor  does  it  require  a  pair  of 
phase-matched  microphones  as  does  method  6.  Additionally,  characteristic 
properties  of  an  absorbing  material  need  not  be  known  as  in  method  4.  Also, 
the  acoustic  properties  of  the  bulk  fiberglass  can  be  determined  using  only 
thin  samples. 

The  theory  developed  for  the  new  method  is  presented  in  Chapter  II  along 
with  a  general  overview  of  the  equations  used  in  the  classical  methods.  The 
experimental  results  for  each  of  the  methods  used  in  comparison  are  presented 
in  Chapter  III  and  comparisons  are  made  and  discussed  in  Chapter  TV. 
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TABLE  1 

ADVANTAGES  AND  DISADVANTAGES  OF  METHODS  USED  TO  DETERMINE  THE  CHARACTERISTIC 

IMPEDANCE  AND  COMPLEX  l^VE  NUMBER 


Method 

1 .  Probe  Tube 
(Scott) 


Advantages 

a)  Simple  theory. 


b)  Simple  procedure, 


Disadvantages 

a)  Thick  sample  needed 
to  get  character- 
istic impedance 
without  a  wave 
model  for  the 
porous  medium. 

b)  Must  probe  sound 
field. 


2 .  Impedance  Tube 
( Beranek ) 


a)  Simple  theory. 


a)  Thick  sample  needed 
to  get  character- 
istic impedance 
without  a  wave 
model  for  porous 
medium. 


b)  Simple  procedure. 

c)  Probing   of  material 
not  reguired. 

d)  Normal  absorption  easy 
to  measure. 


3 .  Resonance 
( Beranek ) 


a)  Probing  not  required. 


a)  Theory  is  somewhat 
difficult. 


b)  Measurement  of  re- 
sonance frequency  and 
half  power  widths  is 
precise. 


b)  Discrete   frequen- 
cies. 


c )  Measures    surface 
impedance . 


4.  Material  Properties 
(Beranek) 


a)  Sound  wave  not  needed. 


d)  Requires  low  losses 
for  well  defined 
resonances . 

a)  Complex  theory. 

b)  Measurement  of 
microscopic  proper- 
ties is  complex. 
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TABLE  1  (continued) 


net hod 


Advantages 


Disadvantages 

c)  Requires  classifi- 
cation of  material 
as  rigid  or  soft. 


5.  Yaniv 


6.  Two  Point  Phase  Diff- 
erence (Chung  and 
Blaser) 


a)  Simple  theory. 

b)  Single  grocedure. 


c)  Probing  of   material 
not  necessary. 

a)  Simple  theory. 


b)  Simple  procedure. 


a)  Time  consuming. 

b)  Quarter  wavelength 
termination  needs 
continuous  adjust- 
ment with  freguen- 
cy. 


a)  Thick 
needed . 


sample 


b)  Microphone  posi- 
tions only  optimum 
for  select  freguen- 
cies. 


c)  Probing   of  material 
not  necessary. 


c)  Reguires  accurate 
microphone  nhase 
calibrations. 


7.  New  Method 
(Baker) 


a)  Probing  not  required. 


b)  Simple  procedure. 


a)  Complex  theory  and 
numerical  methods 
required. 


b)  Discrete 
cies. 


frequen- 


c)  Ability  to  determine 
DC  flow  resistance 
non-invasively. 

d)  Bulk  properties 
obtained  from  thin 
samples. 


c)  Requires  low  losses 
for  well  defined 
resonances . 
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II.   THEORY 

The  theory  behind  the  new  method  is  presented  in  detail  in  this  section 
with  much  of  the  fundamental  concepts  and  understanding  havinq  cone  fron  the 
investigation  of  sound  propagation  in  a  porous  solid  saturated  with 
super- flu id  Helium  [Ref  11] .  Additionally,  the  basic  theory  and  eguations  used 
are  presented  for  each  of  the  classical  methods  used  in  comparison  with  the 
new  technigues. 


A.   THEORY  OF  A  FLUID-SATURATED  RIGID  POROUS  SOLID 

The  linearized  eguation  of  motion  for  the  microscopic  fluid  velocity 


'micro 


in  the  Eulerian  description  is  [Ref.  12,  p.  49], 


Pf 


3v, 


micro 


3t 


+  '  micro*  ^vmicro 


(2.1) 


■*      ■*■  -i    •>  ■*  > 

=  -Vo  +  nV^  •     +  (r  +  —   n)7(V'v  •    ). 
M    '   micro   ^  -t       u    v   micro' 


Using 


V2v  .   =  v(v«v  •   : 

v  micro    K        micro' 


-  V  x  V  x  v  • 

micro 


(2.2) 


we  can  write  Eg.  (2.1)  as 
^micro 


Pf 


at 


+  (vmicro-V)v, 


micro 


=  -Vp  +  nV  x  V  x 


+  U  +  5  n)V(V.vnirr 

(2.3) 


micro   v ^   3  '"  "x  '  micro 


Where  pf  is  the  fluid  mass  density,  p  is  the  fluid  pressure,  n  is  the  fluid 
shear  viscosity  and  5  is  the  fluid  bulk  viscocity.  The  last  term  in  this 
eguation  can  be  ingnored  here  because  its  affect  is  small  compared  to  the 
other  terms  on  the  scale  of  a  pore.   Ignoring  the  last  term,  Eg.  (2.3)  becomes 


Pf 


Dv 


micro 


Dt 


=  -Vp  -  nVx7xv, 


micro  ' 


2.4) 
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f^niicro 
where  — is  the  total  time  deriviative  of  the  microscopic   fluid 

velocity.  From  Eq.  (2.4),  a  phenomenological,  linearized  equation  for  the 
average  (over  the  microscopic  structure)  fluid  velocity,  denoted  <vmicro>'  can 
be  written: 

Peff    dt  =  "V<P>  "  Rflow<Vmicro>  •  <2-5> 

Here  Peff  i-s  ^e   phenomenological  effective  fluid  mass  density  and  Rf]_ow  is 

the  phenomenological  flow  resistance  [Ref .  13,  pp.  252-255] .  Rflow  ^s  defined 

such  that  in  steady  state 

-V<p>  =  Rqc  <^micro>   •  (2.6) 

In  general  (including  oscillatory  flow)  peff  and  Rf]_ow  can  be  operationally 
defined  by 

\   Peff  <Vcro>2    =  ^Cfc^vofSme"57  per  mit    (2*7) 

D     x-v     •> 2       =  R^te  of  dissipation  of    (2.8) 
Rflow  ^  micro  kinetic  energy  density 

For  eJ^  time  dependance,  Eg.  (2.5)  can  be  written,  upon  combining  both 
velocity  terms  on  the  left  hand  side, 

<lirrn>    =   "    V<p>       ,  (2.9) 


micro 


ju  f%ff(w)  +  Rfiow(w) 

where  peff  and  Rf]_ow  are  now  freguency  dependent  phenomenological  properties 
which  describe  the  effect  on  the  fluid  due  to  the  microscopic  geometry  [Ref. 
II,  pp.  334-335].  It  should  be  noted  that  in  Eq.  (2.9)  and  in  the  remainder  of 
this  thesis  complex  quantities  are  denoted  by  a  tilde. 
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Defininq  complex  effective  density,  peff(Gj),  as 

peff(u)  -  peff(u))  + — j  (2.10) 

J  w 

Eq.  (2.9)  nay  be  written  in  a  nore  compact  form  as 

^       -*  _ 

jajpeff(w)  <Vnicro>  =  "7<p>   ■  (2.11) 

Plane  progressive  wave  solutions  to  Eq.  (2.11)  travellinq  in  the  positive 
x-direction  nay  be  souqht  by  lettinq  the  acoustic  pressure  be  qiven  by 
Sp  eJ(tjt  "  kx)  and  lettinq  <vmicro>= vavq  eJ(wt  -  kx),  where  k  =  k  -  ja  with 
k,  a  >  0  such  that  the  enerqy  of  the  wave  exponentially  decays  as  one 
moves  deeper  into  the  porous  solid  (note  in  order  to  avoid  confusion  later, 
the  letter  a  will  be  used  to  denote  the  attenuation  constant,  reserving 
the  letter  a  for  the  tortuosity  in  the  Biot  theory).  Eq.  (2.11)  then 
becomes 

jw  p  C£(w)  v   =  j  k~  £p  .  (7  -ij) 

J       eff     avg  i^.izj 

The  microscopic  equation  of  continuity  for  vnicro  stiH  h°lds  f°r  <Vmicro> 

if  we  consider  a  larqe  enouqh  volume  so  that  irreqularities  on  the  order  of 

the  pore  size  averaqe  out,  and  if  we  consider  the  density  of  the  fluid  in  the 

pores  to  oscillate  about  the  bulk  fluid  density  with  the  passing  of  an 

acoustic  disturbance.  Hence  we  let  p  =  pf  +  5p  el'wt  "  kx)  and   the  equation 

of  continuity  becomes,  to  first  order, 

jufip  -  jpf  k  vavq  =  0   .  (2.13) 

The  relation  between  pressure  and  density  variations  (the  equation  of 
state)  in  a  fluid-saturated  porous  solid  can  be  written 

#=  cf2  3£»  ,  (2.14) 

Sp     f    Y 
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where  Cf  is  the  speed  of  sound  in  the  bulk  fluid,  y(w)  is  the  complex  ratio 
of  specific  heats  (Cp/cv)  as  a  function  of  frequency,  and  y  is  the  adiabatic 
ratio  of  specific  heats. 

Eq's.   (2.12),   (2.13)   and  (2.14)  lead  to  a  wave  equation  with  the  phase 
velocity  determined  as  follows:  writinq  Eq.  (2.12)  (the  equation  of  motion) 
as 

|  Peff^avg  =  -f§  ^  (2.15) 

and  substituting  Eq.  (2.14)  into  it  results  in 

^  Peff(^)  vavg  "  cf  ^ 5P  •  (2-16> 

Recastinq  Eq.  (2.13)  (continuity)  as  pf  vavg  =M  6p  and  dividing  Eq.  (2.16) 
by  it  results  in 

l^r  =  cf  "Y"-  '  Pfeff{S)  (2'17) 

Pf 

Defininq  the  complex  tortuosity  in  the  fluid  saturated  porous  solid  by  [Ref. 
8,  p.  185  and  Ref.  11,  p.  37] 

a(w)  =  -^ ,  (2.18) 

Pf 

the  phase  velocity  can  be  written  as 


<-  (#)2=Scr  ^  .  (2.19) 

VK  ;     a(  oj)    Y 


The  complex  tortuosity  can  be  written  as 


RnrF(u) 

o(u)  =  a  +  -=^ (2.20) 

Jwpf 
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where  a  is  the  intrinsic,  or  infinite  frequency  (inviscid)  tortuosity,  and 
F(  oj)  is  a  dinensionless  correction  function  whose  purpose  is  to  make  the 
frequency  dependence  of  a( w)  cone  out  riqht  [Ref.  8,  p.  185  and  Ref.  11,  p. 
37]. 

The  complex  tortuosity  can  be  written  in  the  form  [Ref.  11,  p.  49] 


where 


aU)  =  a  -  jrFU)   ,  (2.21) 


r  =  ^  (2.22) 

to  pf 

is  a  dimensionless  flow  resistance. 

The  mathematical  form  for  the  dimensionless  correction  function  can  be 
taken  as  [Ref.  11,  p.  57  and  Ref.  14] 

FU)  =   1  +  j(|)2       ,  (2.23) 

where  ~  .._ 

<  =  5(-)"     .  (2.24) 

Here  5  is  the  structural  factor  in  Biot's  universal  correction  function  and  a 
is  the  intrinsic  tortuosity  of  the  microscopic  pore  qeometry  [Ref.  10] . 

The  intrinsic  tortuosity  is  a  weiqhted  path  lenqth  correction  factor  which 
satisfies  a>l  .  It  depends  most  stronqly  on  porosity  P  ,  with  a*l  as 
P+l  (as  is  the  case  for  fiberqlass).  Typical  values  ranqe  from  1  to  4  [Ref. 
11,  p.  227  and  Ref.  15] . 

The  DC  flow  resistance  is  proportional  to  n/(pore  size)2,  so  that 

r  _  n         1       -  .viscous  penetration  depth.  2      ,?  ?c-  ] 
a   aupf  *  (pore  size) l        l        pore  size  * 
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The  quantity  r/ct  determines  the  relative  importance  of  flow  resistance  versus 
fluid  inertia  in  determining  the  microscopic  flow  field. 

The  structural  6  depends  upon  the  averaqe  cross-sectional  shape  of  a 
pore  [Ref.  8,  p.  185],  Experimentally  it  appears  that  5h4  for  many  qranular 
porous  solids  [Ref.  16]  . 

Substituting  Eq.  (2.22),  (2.23)  and  (2.24)  into  Eq.  (2.21)  results  in 

ZiZ)    =a^i    t^k^+i    (|)2  [0)11/21   •            (2-26) 

For  all  calculations  in  this  research  it  will  be  assumed  that  a=l  and  5=4  for 
fiberqlass   [Ref.  11,  pp.  49-59]. 

The  complex  ratio  of  specific  heats,  y(ui)  ,  is  defined  as 

y(gj)  =  Cp/cv  (2.27) 

where  cD   and  cv  are,  respectively,  the  complex,  frequency  dependent  specific 

heats  per  qram  at  constant  pressure  and  at  constant  volume.  For  an  ideal  qas 

~  ~         R 

cD  -  Cy  =  —  ,   where  CI   is   the  molar  mass,   and  so 

YU)    =  1  +  pjf-     •  (2.28) 

The  propaqation  of  sound  in  the  porous  medium  will  vary  from  isothermal 
at  low  frequencies  when  the  thermal  penetration  depth  is  larqe  compared  to 
a  pore  size,  hence  cv  +  °°  ,  to  adiabatic  at  hiqh  frequencies  when  the  thermal 
penetration  depth  is  small  compared  to  a  pore  size,  and  cv  =  bulk  fluid 
value.  Since  the  transport  of  momentum  and  enerqy  in  a  qas  both  arise  from 
molecular  collisions,  it  will  be  assumed  that  thefunctional  form  (y  -  1)  has 
the  same  form  as  the  complex  tortuosity  a(  a))  ,  only  with  the  thermal 
penetration  depth  5^  in  place  of  the  viscous  penetration  depth  5V^SC. 
Hence  we  take  [S.  Baker,  private  communication] 

7(5  =1  +  ^1  (2.29) 
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where  r  and  F(  a>)  are  given  by  Eq's   (2.22)   and   (2.23),   respectively.  The 
viscous  penetration  depth,  5visc  '  and  tne  thermal  penetration  depth  5th  are 

defined  by 


Svisc  =  J—  «th  =  \    ^h~  (2-3°) 

where  n  is  the  fluid  shear  viscosity,  pf  is  the  bulk  fluid  density,  <  is 
the  thermal  conductivity  and  cp  is  the  specific  heat  per  gram  at  constant 
pressure. 

4 

The  factor  — can  be  recast  as 

'Vise 

^  •  -  *    •    i  (2.31) 


2 

6 


vise 


where 


and 


v  =  n/pf  =  kinematic  viscosity  (2.32) 


X  =  k/p  cp  h  thermal  diffusivity  (2.33) 

and  Pr  =  Prandtl  number  which  has  a  value  of  0.733  for  air  at  20°C.  Hence  the 
complex  ratio  of  specific  heats  can  be  rewritten  as 

7(~)  =1+I^Uf^   .  (2.34) 

Pr 

In  the  high  freguency  limit,  the  sguare  of  the  phase  velocity  in  the 

porous  material  (Eg.  2.19)  goes  to  the  fluid  inertia  dominated  limit  in  which 

sound  propagates  at  the  bulk  fluid  speed  reduced  by  the  sguare  root  of  the 

tortuosity.   Conversely,   in  the  low   freguency  limit,   flow   resistance 
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dominates,  and  pressure  oscillations  diffuse  throuqh  the  porous  material.  The 

crossover  from  inertia  dominance  to  flow  resistance  dominance  occurs  when  the 

viscous  penetration  depth  is  approximately  equal  to  a  characteristic  pore 

size.  Figure  2.1  summarizes  the  limiting  conditions. 

1 .  Development  of  Dispersion  Relation  in  a  Fluid-Filled  Cavity  Partially 
Filled  With  a  Porous  Solid 

Consider  a  rectangular  cavity  which  is  partially  filled  with  a  fluid- 
saturated  porous  solid.  The  region  occupied  by  bulk  fluid  is  denoted  as 
region  1;  the  region  occupied  by  the  fluid-saturated  porous  solid  is 
identified  as  region  2.  The  porous  solid  is  assumed  homogeneous  and 
isotropic.  Figure  2.2  depicts  the  geometry  used  to  develop  the  dispersion 
relation  for  standing  waves  in  the  cavity.  It  will  be  assumed  that  all  cavity 
cavity  walls  are  perfectly  rigid  so  that  n«u  =  0.  Then  n*Vp  =  0  at  these 
boundaries  and  the  velocity  potential  in  regions  1  and  2  can  be  written  as 
follows: 

In  region  1  (bulk  fluid), 


where 


and 


l  =   <J10  e-1"     cos[kxx]    costk^U-d!)]  (2.35) 


k     =  j1  ;        n  =  0,    1,    2,    .    .    .  ,                                     (2.36) 

Z2  =    (kx2  +  klz2)    cf2  (2.37) 

cf2  =  Bf/pf        ,  (2.38) 

Bf   =  bulk  modulus  of  fluid 

<j>10     =  complex     amplitude     of     the     velocity  potential     in     region     1. 
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High  and  Low  Frequency  Limits  on  Phase  Velocity 
in  a  Porous  Material 


High  Frequency  Limit 


! 


Low  Frequency  Limit 


"SI 


=5^ 


Plug  Flow 


Laminar  Flow 


S(  td)   +  Ooo 


a(oj) 


pf3 


Y(  u>)  +  y 


YU)  *  1 


Cf 


•  Cf  "Pf 
Y  KDC 


+  Fluid  Inertia  Dominates 
■>  Sound  Propagates 


♦  Flow  Resistance  Dcminates 

■*■  Sound  does  not  Propagate 
(Pure  Diffusion) 


Cross-over  when: 
^visc  a  characteristic  pore  size 


Figure  2.1  High  and  Low  Frequency  Limits  on  Phase  Velocity  in  a  Porous 
Material 
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REGION  JL 
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•                                        •                          » 

.    .          ... 

•                                                             •                • 
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1-  -«*t 
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Figure  2.2  Rectangular  Cavity  Geometry  for  New  Method, 
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The   x-conponent     of    the  wave   number    is     a  pure    real   number  since      the   boundary 
conditions  of    rigid  walls   normal      to   the   x-direction  at   x     =0     and   x  =     L   force 


real  and  discrete  values   for  kx   . 


In   reqion   1    the  particle  velocity  can   be  written  as 

Vj   =    7^  (2.39) 


and   the  pressure  as 


where 


Spl  ~  -jw  Pf    <fl   •  (2.40) 

In   reqion  2    (porous  material), 

*2  =   *20  e  ^  cos[k  x]    cos[k2Z(z+d2)]  (2.41) 

A. 


kx  =  ~     .        n  =  o,    1,    2,    .    .    .    ,  (2.42) 

2         „      2    _    2 
^2   =    (k         +    k2z    )Cz      t  (2.43) 

X. 

~    2        cf  yU)  ,.     ... 

c    2  =  ^T^V       >  (2.44) 

L         a(  w)         y 


and 


4>2o  =  complex  amplitude  of   the  velocity  potential    in  reqion  2. 
The     particle  velocity   in  reqion  2   is 

v2  =    1%2  (2.45) 

and   the  pressure    is,    from  Eq.    (2.11), 

6p2  -  -ju)a(co)    pf    <j>2   •  (2.46) 
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With  the  velocity  potentials  defined  in  both  regions  1  and  2  all  that 
is  left  is  to  match  conditions  at  the  boundary  between  the  bulk  fluid  and 
porous  material  (i.e.,  at  z  =  0)  to  produce  the  desired  dispersion  relation. 
The  two  boundary  conditions  that  must  be  satisfied  at  z  =  0  are: 
1.  continuity  of  pressure  across  the  boundary: 


(J>1(z=o)  =  a(u>)  <j>2(z=o) 


(2.47) 


and 


2.  continuity  of  the  normal  component  of  volume  velocity: 


3*1 


=  P 


H2 


z=o 


z=o 


2.48) 


where     P     is   the     porosity  of     the     material.     Applyinq     boundary  condition     1 
results   in 


<j>10  costkizdj   =  ct(u))    4>20  cos[k2zd2] 


Applying  boundary  condition  2  results  in 


(2.49) 


<j>10  klz  sin    [kizdj      =     -J20  p  k2z  sin[k2zd2]      •  (2.50) 

Dividing  Eg.    (2.50)   by  Eg.    (2.49)    results   in  the  dispersion  relation 


where 


klz  tan     [kizdj   +     ~(~>    k2z  tan    [^^2!    =  °     1  (2.51) 


1Z 


UJ  , 


1/2 


(2.52) 
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*2Z 

and 


~2  2 

c2   aU)  ~WZ)~   kx 

f 


1/2 

(2.53 


kx  =  ^  ;   n  =  0,  1,  2,  .  .  .   .  (2.54 


In  Eq.  (2.53)  the  forms  of  a(w)  and  y( u)  have  already  been  defined  in  Eq. 
(2.26)  and  Eq.  (2.34),  respectively.  Eq.  (2.51)  is  a  transcendental  equation 
in  jj  whose  roots  determine  the  natural  frequencies  and  half  power  widths  or, 
equivalently,  the  quality  factors,  of  the  allowed  "standinq  wave"  modes  of 
this  system.  From  these  roots  can  be  determined  the  complex  wave  numbers 
k-^z  and  k   ar>d  the  result inq  pressure  distribution. 

Two  approaches  to  the  use  of  the  dispersion  relation  become  evident. 
The  first  is  the  "forward"  problem,  where  the  roots  of  Eq.  (2.51)  for  oj  may 
be  reqarded  as  functions  of  Rqc  and  a  ,  and  measured  resonance  frequencies 
and  Q's  may  be  compared  with  predictions;  the  second  in  the  "inverse" 
problem,  where  R^  and  a  are  reqarded  as  functions  of  the  roots  of  Eq. 
(2.51)  for  to  and  are  found  by  variation  such  that  the  best  fit  is  obtained 
between  the  measured  and  calculated  spectra.  The  aim  of  this  work  is  to 
verify  the  technique  to  solve  the  inverse  problem  by  oomparinq  the  value  of 
Rqc  extracted  from  resonance  data  with  that  actually  measured.  Appendix 
A  discusses  the  numerical  technique  to  determine  the  roots  of  the  dispersion 
relation  and  the  alqorithm  used  to  determine  the  best  fit  value  of  Rqc* 

Once  the  "best  fit"  value  of  R^c  has  been  determined,  then,  since 
the  roots  of  the  dispersion  relation  in  terms  of  u  are  known  for  this  value 
of  Rqc  /  the  complex  wave  number  is  found  from 

k  =  |  =  k  -  ja  (2.56) 

c 
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with  c  found  from  Eq's.   (2.19),   (2.26)   and   (2.34).  The  characteristic 
impedance  can  be  calculated  from 

2char  =  feff<S>  ?  <2-57) 

where  ^eff     ^s  Qiven  bY  EQ*   (2.10).   Noramlizing  the  characteristic 
impedance  to  that  of  the  bulk  fluid  (/°oc)  results  in 

2?char  (normalized)  =  Zcnar//>0c   .  (2.58) 

The  complex  pressure  field  for  both  regions  may  be  calculated  once 

the  roots  of  Eq.  (2.51)  are  found.  The  pressure  distribution  in  region  1  can 
be  written  in  terms  of  the  velocity  potential  in  reqion  1  as 

£p,  =  -fift  ?.  '  <2-59> 

where  <f>  i   is  given  by  Eq.  (2.35).  The  pressure  distribution  in  region  2  can 
be  written  in  terms  of  the  velocity  potential  in  region  2  as 

(2.60) 


.  *»  "*  ..«•' , 


where  ^2  ^s  Qiven  by  EQ»  (2.41). 

The  velocity  potentials  <p  \  and  02  are  related  by  the  boundary 
condition  that  pressure  be  continuous  across  the  interface  between  the  bulk 
fluid  and  the  porous  material,  Eq.  (2.47),  which  requires  that 

0IO  cos[klZd,]  =  «(uj)02O  cos[k2zd2]   .  (2.61) 

For  simplicity,  choose  0  20  to  depend  on  010  '  resultinq  in 

2  0io   cos  k,zd, 

^rt  =^7^7  *» h     .  (2.62) 

20   oc(u>)  cos  k2Z  djj 

The  pressure  distribution  in  each  reqion  may  be  written  as  follows: 


a. 


P. 


-jco  y»£  cos[k  x]    cos[k,  z(z-d( )] 

2\ 
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?l0        eJWt,  (2.63) 


and  for  region  2, 


cos  klzd 
-j2  pf  cos[kxx]    cqs  ^     d       cos[k2z(z+d2)] 


<J>10  eJ     (2.64 


With  Eq.  (2.63)  and  Eq.  (2.64),  and  the  values  of  k]_z  and  k2z  determined 
from  the  roots  of  the  dispersion  relation  of  Eq.  (2.51),  the  pressure 
distribution  throughout  the  cavity  may  be  calculated.  See  Appendix  B  for  a 
more  detailed  discussion  of  the  pressure  calculation. 

B.   CURRENT  METHODS  USING  TRANSMISSION  LINE  THEORY 

Three  of  the  five  classical  methods  for  determininq  the  characteristic 
impedance  and  complex  wave  number  of  a  porous  acoustic  material  are  based  on 
transmission  line  theory.  These  methods  involve  measurement  of  various 
properties  of  the  stand inq  wave  and  data  are  taken  at  various  points  somewhat 
removed  from  the  sample  surface  [Ref.  4,  pp.  317-321].  The  three  methods 
considered  are: 

—  Impedance  Tube  Method 

—  Resonance  Analysis  Method 

—  Yaniv's  Method 

Only  Yaniv's  method  allows  direct  determination  of  the  material's 
characteristic  impedance  and  complex  wave  number.  The  impedance  tube  method 
allows  direct  measurement  of  the  resistive  and  reactive  components  of  the 
surface  impedance.  The  complex  wave  number  can  only  be  calculated  if 
reflections  from  the  far  side  of  the  porous  material  can  be  neqlected  and  one 
has  some  knowledqe  of  the  complex  effective  mass  density  of  the  fluid  in  the 
pores.  As  developed  by  Beranek,  the  resonance  anlaysis  method  measures  the 
same  quantities  as  the  impedance  tube  method,  with  the  same  limitations.  It  is 
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discussed  primarily  for  its  similarity  to  the  new  technique  in  that  the 
parameters  of  the  resonance  curve  are  measured. 
1 .   Impedance  Tube  Method 

In  this  method  the  non-dissipative  acoustic  transmission  line  is 
considered  a  smooth  riqid-walled  tube  whose  lateral  dimensions  are  small 
compared  to  a  wavelength.  The  losses  at  the  wall  and  in  the  bulk  volume  are 
considered  small  compared  to  the  losses  associated  with  the  absorbing  material 
under  test.  Assuming  that  the  source  of  sound  is  at  one  end  of  the  tube  (x=0) 
and  the  sample  under  test  is  located  at  the  opposite  end  (x=L),  the  surface 
impedance  can  be  determined.  Figure  2.3  depicts  the  geometry  of  the  impedance 
tube. 

The  impedance  tube  that  has  a  constant  cross  section  S  and  the 
source  vibrates  harmonically  at  a  freguency  sufficiently  low  that  only  plane 
waves  can  propagate.  A  standing  wave  pattern  is  set  up  by  the  right  and  left 
travelling  waves  and  can  be  expressed  mathematically  as 

~=  ~eJ(Ut  +  k(L-x)]  +geJ[ut-k(L-x)]  (2_65) 

where  p  is  the  acoustic  pressure,  A  and  B  are  determined  from  the  boundary 
conditions  at  x=0  and  x=L. 

It  has  been  shown  [Ref.  17,  pp.  200-206]  that  the  surface  impedance 
of  the  termination  may  be  written  in  the  form, 

2L    1  +  B/A  ej9 


P°C   1  -  B/A  eJ 


(2.66) 


where  B/A  is  the  ratio  of  the  relected  wave  amplitude  to  the  incident  wave 
amplitude  and  0  is  the  phase  difference  between  the  relected  and  incident 
wave.  Thus  given  B/A  and  0  ,  the  surface  impedance  of  the  termination  can  be 
calculated. 

31 


Solving  for  the  amplitude  |p|  of  the  standing  wave  it  can  be  seen  that 

the  amplitude  at  pressure  antinode  is  A+B,  and  at  a  pressude  node  is  A-B.  The 

ratio  of  the  pressure  amplitude  at  an  antinode  to  that  at  a  node  is  the 

standing  wave  ratio,  or  SWR: 

™m   A+B     Panti  ,_  __. 

SWR  =  T-r-  =         ,  (2.67) 

A"B     Pnode 

which  can  be  rearranged  to  give  the  ratio  B/A: 

B  _  SWR-1 

A  "  SWR+1  *  (2'68) 

Thus  measurement  of  the  SWR  by  probing  the  sound  field  in  the  impedance  tube 

with  a  small  microphone  yields  the  value  of  B/A  [Ref .  17,  p.  206] . 
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Figure  2.3  Geometry  for  the  Impedance  Tube 

The  phase  angle  9  can  be  evaluated  from  determining  the  distance  to 
the  first  node  from  the  end  at  x=L.  These  nodes  are  located  at 

k(L-xn)  -  6/2  =  (n-l/2)it  ;  n  =  1,  2,  .  .  .   ,  (2.69) 

so  that  the  phase  angle  for  the  first  node  is 

Ql   =  2k(L-x1)  -  tt  .  (2.70) 
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With  the  SWR  and  9]_  known,  the  resistive  and  reactive  portions  of  the 
impedance  can  be  determined  graphically  using  a  "Smith"  chart.  To  determine 
the  characteristic  impedance  of  the  porous  material,  the  sample  thickness  is 
increased  such  that  the  reflected  wave  from  the  rigid  termination  is 
insignificant  compared  to  the  incident  wave.  From  the  definition  of  the 
material's  characteristic  impedance,  the  complex  wave  number  can  be  found 
assuming  that  we  know  how  the  fluid  reacts  in  the  pores  of  the  material.  The 
characteristic  impedance  may  be  written  as 

zchar  =  Peff  c  (2'71) 

where  peff  is  the  complex  effective  fluid  density  in  the  pores  and  c  is  the 

phase  velocity  of  the  acoustic  distrurbance  in  the  pores.   The  phase  speed 

c  in  a  fluid,  can  expressed  as 

c2  =  Bf/peff  (2.72) 

where   EL-  is  the  complex  bulk  modulus  of  the  fluid  in  the  pores.   Therefore 

Eg.  (2.71)  becomes,  upon  substituting  Eg.  (2.72)  for  Peff' 

Zchar  =  Bf/c  (2.73) 

and  since  c  =  =-,  the  eguation  that  describes  the  complex  wave  number  K  can 
be  written  as 

k  =  $£  Zchar  •  <2-74) 

If  is  is  assumed  that  the  material  pore  size  is  small  compared  to  the  viscous 
and  thermal  penetration  depths  (in  air  these  two  parameters  are  approximately 
egual)  in  the  freguency  range  of  interest,  we  can  write 

Bf  =  Bp  (2.75) 

where  Bp  is  the  isothermal  bulk  modulus  of  the  fluid  in  the  pores. 
Additionally  for  an  ideal  gas,  the  isothermal  bulk  modulus  can  be  written  as 
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Bj-  =  1/ttj.  (2.76) 

where   <>p  is  the  isothermal  compressibility  of  the  gas,  and  is  equal  to  1/P, 

where  P  is  the  pressure  in  the  pores  [Ref.  18,  p.  47] .  Hence  Eq.  (2.74)  can 

be  written,  assuming  isothermal  conditions  exist,  in  the  form, 

k  =  uitfp  Zchar  (2.77) 

where  k  =  k  -  ja,  for  a  wave  travelling  in  the  positive  x-direction  in  the 
porous  solid.  Eg.  (2.77)  is  presented  to  show  the  relationship  between  the 
complex  wave  number  and  the  characteristic  impedance  assuming  isothermal 
conditions  in  the  pore  volume.  This  equation  was  used  only  to  check  Beranek's 
Material  Properties  Method  in  the  isothermal  limit.  Only  the  values  for 
^char  extracted  by  the  impedance  tube  method  are  reported  here. 
2 .  Resonance  Analysis  Method 

The  resonance  analysis  method  is  used  to  determine  the  surface 
impedance  of  a  porous  acoustic  material.  This  method  reguires  low  loss  within 
the  material  to  be  successful.  A  resonance  curve  may  be  traced  out  by  varying 
either  the  tube  length  or  the  freguency  in  this  method.  For  the  present 
research,  the  impedance  tube  discussed  previously  was  used,  with  freguency 
swept  to  trace  out  a  resonance. 

Beranek  presents  the  eguation  needed  to  determine  the  impedance  of  a 
sample  under  test  and  is  rewritten  here  only  for  discussion  [Ref.  4,  pp.  329- 
334]  . 

The  normalized  surface  impedance  of  the  sample  is  written, 

l?u£f  m   y/j  +  k^  a  cQth  [(   _  j^ji/c]         (2.78) 
Poc 

where  v  =  |Zsurf |/pQc  ,  $  being  the  phase  angle  of   |Zsurf |  and  Zsurf 
is  the  impedance  at  the  surface  of  the  sample  which  is  to  be  determined.   At 
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the  resonance  peak,  the  peak  frequency  at-,  is  approximately  equal  to  the 
normal  angular  frequency  of  lonqitudinal  vibration,  u^  ,  or 

ajp  =  a^   .  (2.79) 

At  the  half  power  frequencies  (-3  dB  down  from  the  peak)  the  half-power  half 
width  (HPHW)  can  be  written  as 

HPHW  =  ^  "  'H   =  ^-*I  a  kx   ,  (2.80) 

Z         Up 

where  k^  =  k^  +  kna,  and  kna  is  the  spatial  damping  constant  without  the 
porous  material  present  (i.e.,  it  takes  into  account  wall  and  driver  losses). 
Therefore  km/  the  spatial  damping  factor  due  to  the  sample  under  test  can  be 

written  as 

k  =  k i  _  k  (2.81 ) 

m    l         na  u,0i' 

Assuming  the  ratio  k^o^  <<  1/  the  surface  impedance  of  the  sampole  can  be 
determined  once  u^,  w^r  ojp,  I,   and  c  are  measured. 
3.  Yaniv's  Method 

Equations  have  been  developed  in  transmission  line  theory  for 
determination  of  the  input  impedance  of  an  electrical  transmission  line  in 
terms  of  the  line  characteristics:  propagation  constant,  characteristic 
impedance  and  length  of  line  [Ref .  5] .  When  the  line  is  terminated  in  an 
open  or  short  circuit,  the  input  impedance  of  the  line  may  be  expressed  in 
terms  of  the  line  parameters  with  the  following  equations: 

Zos  =  Zchar  coth  [yL]  (2.82) 

and 

Zss  =  Zchar  tanh  [yL]  (2.83) 

where  Zos  is  the  input  impedance  to  a  line  terminated  in  an  open  circuit 
(Figure  2.4),   Zss  is  the  input  impedance  to  a  line  terminated  in  a  short 
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circuit  (Figure  2.5),  Zchar  is  the  characteristic  impedance  of  the  line, 
Y  =  jk  is  the  complex  propagation  constant  of  the  line  and  L  is  the  length 
of  the  line. 

Fran  Eg.   (2.82)  and  Eg.   (2.83)  we  can  solve  for  the  material's 
characteristic  impedance: 

zchar  "  tz~ss  W  1/2  '  (2.84) 

and  for  the  complex  propagation  constant 

tanh[YL]  =  [2ss/2os\l/1     .  (2.85) 

The  high- impedance  termination  for  an  acoustical  line  may  be  realized  by 

termination  with  a  rigid  wall  and  a  low- impedance  termination  for  an 

acoustical  line  may  be  realized  by  terminating  the  line  with  a  column  of  fluid 

egual  to  a  guarter  wavelength,  itself  terminated  in  a  rigid  wall  [Ref .  5] . 
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Figure  2.4  Geometry  for  Yaniv's  Open  Circuit  Determination 
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Figure  2.5  Geometry  for  Yaniv's  Short  Circuit  Determination 
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By  measuring  the  standinq  wave  inside  the  tube  alternatively 
terminated  with  an  open  circuit  and  then  a  short  circuit,  both  Zos  and  Zgg 
may  be  determined.  With  these  two  values,  Eq.  (2.84)  and  Eq.  (2.85)  may  be 
solved  for  the  desired  complex  propagation  constant  and  characteristic 
impedance . 

C.   PROBE  TUBE  METHOD 

This  method  was  first  presented  by  Scott  (1946)  [Ref.  2]  to  allow  the 
determination  of  the  complec  wave  number  from  the  characteristics  of  an 
acoustic  disturbance  propagating  through  a  fluid-saturated  porous  solid  by 
probing  the  sound  field  existing  in  the  material  and  extracting  the  wave 
number  and  absorption  coefficient. 

To  determine  the  complex  wave  number  in  the  porous  medium,  which  describes 
the  spatial  variation  of  the  acoustic  disturbance  for  a  uniform,  progressive 
plane  wave,  Scott  writes  the  velocity  potential  <j>  of  the  average  motion  in  the 
pores  as, 

J  =  x  e-j£x  ejwt  (2.86) 

where  k  =  k  -  ja  is  the  complex  wave  number  and  x  represents  distance  in  the 
direction  of  propagation  of  the  wave.  Eq.  (2.86)  may  be  rearranged  such 
that, 

1     -       «"ax  J(wt  ~  kx>  (2.87) 

<P   _  9oe     e 

where  a  is  the  absorpotion  coefficient  and  k  =  2tt/X  is  the  wave  number. 
Experimentally,  one  measures  the  acoustic  pressure  amplitude  as  a  function  of 
distance  in  the  material,  along  with  the  wave  length. 

To  determine  the  characteristic  impedance  of  the  porous  material,  Scott 
writes, 
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~        ^ef  f  (  ^  w 

Zchar  = T—  (2'88 

where   peff(w)   is  the  complex  frequency  dependant  effective  density  of  the 

fluid  in  the  pores.  Determination  of  Z^^  assumes  of  course,  that  the 

effective  fluid  density  is  known,  for  example,  from  Biot  theory.   Eq.  (2.88) 

was  used  onlv  to  check  that  the  measured  values  of  k  were  consistent  with  the 

values  of  Zchar  from  the  other  methods.  Only  the  values  of  k  measured  with 

Scott's  method  are  reported  here. 

D.   ?1ATERIAL  PROPERTIES  fETHOD 

This  method  requires  determination  of  five  properties  of  the  material  to 
calculate  the  complex  propaqation  constant  and  characteristic  impedance  of  a 
material.  Beranek  [Ref.  4,  Ch.  19]  discusses  this  method  in  detail  and  only 
that  portion  of  the  theory  applicable  to  the  comparison  with  the  new  technique 
will  be  presented. 

Beranek  has  shown  that  both  the  propaqation  constant  and  characteriustic 
impedance  are  derivable  with  the  aid  of  charts  from  the  basic  physical 
properties  of  the  material  [Ref.  4,  p.  839] . 
These  properties  are: 

(a)  the  specific  flow  resistance 

(b)  the  porosity 

(c)  the  structure  factor 

(d)  the  volume  coefficient  of  elasticity  of  the  air  in  the  interstices 

(e)  the  volime  coefficient  of  elasticity  of  the  skeleton  of  the  material. 
First,  one  must  determine  whether  the  material  under  investiqation  can  be 
considered  a  soft  blanket  or  a  riqid/dense  blanket.  This  determination  is 
made  by  considerinq  the  ratio  of  the  volume  coefficient  of  elasticity  of  the 
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air  in  the  interstices  to  that  of  the  skeleton.   If  this  ratio  is  less  than 

-rv  the  material  is  considered  to  be  a  rigid/dense  blanket   [Ref .  4,  p.  840]  . 

The  fiberglass  material  chosen  for  the  rigid/dense  blanket  criterion  for 

Beranek  provides  an  eguation  for  the  propagation  constant  b  =  jk,  of  such  a 
material,  namely 


~   .-   .    /PQkstY   /.   R^ 

b  =  Jk  =  Du  /-Y—  ^/l  ~  ^j^J   ,  (2.89) 

where 

b   =  complex  propagation  constant  (same  as  Yaniv's  y) , 

k   =  complex  wave  number, 

kst  =  structure  factor  (same  as  Biot's  tortuosity,  a), 

K   =  volume  coefficient  of  elasticity  of  air  in  the  interstices   (i.e., 
the  bulk  modulus,  and  can  be  isothermal,  adiabatic  or  complex), 

r-^  =  dynamic  specific  resistance  per  unit  thickness  (approximately  the  DC 
flow  resistance), 

Pq  =  denisity  of  air,  and 

Y   =  porosity  (same  as  P  in  Biot's  theory). 
By  the  use  of  Figure  2.6  and  Figure  2.7  the  magnitude  and  phase  of  b  can  be 
determined.  The  characteristic  impedance  may  then  be  determined  from 

Z  .    =  -j  ^r  b   .  (2.90) 

char    J  gjY 

This  method  suffers  from  the  fact  that  not  all  of  the  material  properties  can 

be  experimentally  determined  and  thus  reasonable  assumptions  for  their  values 

must  be  made   [Ref.  4,  p.  840] . 
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Figure  2.6  Chart  for  Determining  the  Magnitude  of  the  Propagation  Constant  b 
of  Rigid  Tiles  [Figure  19.3  of  Ref.  4,  p.  843] 
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Figure  2.7  Chart  for  Determining  the  Phase  Angle  of  the  Propagation  Constant  b 
of  Rigid  Tiles  [Figure  19.4  of  Ref.  4,  p.  844] 
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III.   EXPERIMENT  AND  RESULTS 

The  objective  of  this  study  was  to  develop  a  new  technique  to  measure  both 
the  complex  wave  number  and  characteristic  impedance  of  a  fluid-saturated 
porous  solid.  As  with  any  new  technique,  comparison  with  classical  methods  is 
a  qood  basis  on  which  to  judqe  one's  theory.  This  section  will  first  present 
the  experimental  results  of  the  five  classical  techniques  discussed  previously, 
and  then  the  new  technique.  A  comparison  of  measured  data  for  all  methods  will 
be  presented  in  Chapter  TV. 

A.   PROPERTIES  OF  THE  POROUS  SOLID 

The  porous  solid  used  in  all  experiments  was  a  semi-riqid  type  of 
fiberqlass,  much  like  Owens-Corning  Tupe  705  fiberglass,  used  in  the  insulation 
of  refrigeration  type  equipment.  This  material  has  sufficient  rigidity  to 
allow  it  to  be  cut  from  sheets  into  a  variety  of  shapes  and  to  allow  it  to  be 
inserted  into  tubes  without  excessive  compression.  This  property  was  important 
in  the  flow  resistance  and  impedance  tube  measurements  where  a  snug  fit  was 
required. 

Discs  with  nominal  dimensions  3-29/32  inches  in  diameter  and  1-29/32 
inches  thick  were  cut  from  three  sheets  of  this  material,  weighed,  and  measured 
to  determine  the  material's  volume.  From  two  sets  of  measurements,  the  average 
material  density  was  found  to  be  43.0  ±  5.0  kg/m3.  For  comparison, 
Owens-Corning  Type  705  fiberglass  is  reported  to  have  a  denisty  of  96  kg/m3 
[Ref .  5] ,  while  rock-wool,  also  known  commercially  as  Stillite,  is  reported  to 
have  a  density  of  80  kg/m3  [Ref.  2]  . 

The  porosity  of  the  fiberglass  material  was  measured  by  immersing  a  1-inch 

cube  in  a  graduated  beaker  of  water  and  measure ing  the  volume  of  displaced 

water.   The  fiberglass  material  itself  was  assumed  to  be  non- absorb ing.   The 

resulting  porosity  was  found  to  lie  between  0.99  and  1.0. 
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1 .   Flow  Resistance  Measurements 

Next,  measurement  was  made  of  the  specific  DC  flow  resistance  per  unit 
lenqth  of  sample.  For  most  homogeneous,  isotropic  materials  this  flow 
resistance  may  be  measured  with  sufficient  accuracy  by  applying  a  known  steady 
air  flow  across  a  sample  and  measuring  the  resulting  pressure  differential. 
The  actual  method  used  is  similiar  to  that  presented  by  Beranek  [Ref.  4,  p. 
850]  . 

The  quantities  measured  were:   (Ref.  4,  p.  847] 
Ap  =  pressure  differential  between  the  two  faces  of  the  sample 
v  =  velocity  of  the  linear  flow  of  air  through  the  sample  due  to  Ap 
The  pressure  differential  was  measured  with  a  standard  water  manometer  (±0.1 
cm  accuracy)  and  the  velocity  of  linear  flow  was  determined  buy  measuring  the 
time  reguired  to  collect  a  sufficient  volume  of  water  (which  drew  air  through 
the  sample  by  a  siphon)   in  a  graduated  cylinder.  The  experimental  apparatus 
used  is  shown  in  Figure  3.1.  All  tube  connections  were  sealed  with  RTV  Silicon 
Rubber  to  prevent  leaks. 

The  velocity  of  linear  flow  produced  by  Ap  is, 

v-jk  -  (3.1) 

where 

V  =  volume  of  water  collected  in  time  t,  and 

A  =  cross-sectional  area  of  sample. 

Then  the  DC  flow  resistance  R^  is  given  by 

Rnr  =  Ag  a  ApAt  „  pgAhAt  (3  2) 

^   vT    VT      VT 

where   p  is  the  density  of  water,  g  is  the  acceleration  of  gravity,   Ah  is 

the  difference  in  water  height  in  the  two  arms  of  the  manometer,  and  T  is  the 

thickness  of  the  sample.  The  units  of  Rqtj  are  in  Rayl/unit  thickness,  which 

can  by  written  in  CGS  or  MKS  units  gm/cm2, sec/cm  or  kg/m2,sec/m,  respectively. 
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Figure  3.1     dc  Flow  Resistance  Measurement  Apparatus 
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In  the  first  flow  resistance  experiment  seven  discs  were  placed  into 
the  larqe  sample  holder,  each  disc  having  a  diameter  of  9.85  cm  and  a  thickness 
of  4.76  cm.  The  samples  were  made  as  close-fitting  as  possible  but, 
unfortunately,  the  die  used  to  cut  these  discs  had  a  diamter  approximately  0.1 
cm  greater  than  that  of  the  sample  holder.  While  providing  a  tight  fit  along 
the  tune  wall,  this  led  to  a  compression  in  the  material's  thickness  by 
approximately  10%  as  they  were  pushed  into  the  holder.  At  that  time  it  was 
thought  that  the  effect  of  this  compression  would  not  be  significant  enough  to 
require  re-machining  the  cuting  die.  Later  measurements  of  flow  resistance 
using  a  different  holder,  however,  gave  values  which  were  a  factor  of  2  smaller 
than  that  for  flow  normal  to  the  fiberglass  surface. 

Measurements  were  conducted  over  a  limited  range  of  air  flow  velocity 
by  adjusting  the  3/4  inch  Hoke  valve  on  the  water  side  of  the  apparatus  to  give 
the  desired  pressure  differential.  When  steady  state  conditions  were  reached, 
the  volume  flow  was  then  plotted  as  a  function  of  volume  flow  rate.  This  is 
shown  in  Figure  3.2.  Within  experimental  error,  the  variation  in  Ap  with  V  was 
found  to  be  linear  up  to  the  highest  flow  rate  measured,  which  was  many  orders 
of  magnitude  greater  than  expected  for  a  sound  wave.  Since  the  sound  pressures 
appearing  at  the  surface  of  a  sample  are  seldom  greater  than  1000  dynes/cm2  in 
practice  (100  Pa  or  0.5  inches  water)  it  can  be  seen  that  even  at  higher 
pressure  dif ferenitals  the  linear  relationship  betwen  Ap  and  V  is  maintained. 
The  DC  flow  resistance  had  a  value  of  59,100  ±  4,200  MKS  Rayl/m  using  the 
uncompressed  thickness.  Table  2  lists  the  experimental  data  and 
Table  3  lists  the  DC  flow  resistance  results  using  the  compressed  and 
uncompressed  thickness. 

It  was  discovered  later  (while  attempting  to  reconcile  these 
measurements  with  the  DC  flow  resistance  required  to  fit  the  resonance  data) 
that  the  DC  flow  resistance  of  our  fiberglass  samples  is  guite  anisotropic  (the 
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DC  flow  resistance  required  to  fit  the  [100]  modes  resonances,  which  has  its 
particle  velocity  parallel  to  the  air-fiberglass  surface,  was  between  10,000  to 
20,000  MKS  Rayl/m) . 


7H.fr 


Press.  Difference  vs  Volume  Flow 

(VOLUME  FLOW  NORMAL  TO  SURFRCE) 


33.1 


VOLUME  FLOW  RATE 
( cc/sac ) 


Figure  3.2   Ap  vs  V  in  the  Large  Sample  Holder 
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TABLE   2 


LARGE   SAMPLE   HOLDER   DC   FLOW  RESISTANCE   MEASUREMENT   DATA 


(Flow  Normal   to  Surface) 


Uncompressed  thickness  =33.34    ±  0.1   cm 
Compressed  thickness  =29.80    ±  0.1  cm 


Sample  holder  diameter  =  9.84    ±0.1  cm 
Room  temperature  =  23.3  deg.  Centriqrade 


Ah 

Vol. 

t 

V 

%:* 

(cc) 

(cc) 

(sec) 

(cc/sec) 

(MKS   Rayl/m) 

6.4    ±  0.1 

2038    ±  4 

7.3    ±  0.1 

278.0    ±  0.5 

57.5xl03    ±  700 

6.4 

2102 

7.6 

277.7 

57.6 

6.3 

2074 

7.4 

278.8 

56.5 

5.6 

2016 

8.4 

240.6 

58.2 

5.5 

2056 

8.7 

236.6 

58.1 

5.3 

2008 

8.9 

225.1 

58.8 

4.5 

2046 

10.3 

198.6 

56.6 

4.5 

2122 

10.6 

201.1 

55.9 

4.3 

2108 

11.4 

184.9 

60.8 

3.5 

2128 

14.1 

151.5 

57.7 

3.5 

2063 

14.3 

144.7 

60.4 

3.5 

2074 

13.5 

153.6 

56.9 

2.5 

2077 

18.1 

114.6 

54.5 

2.5 

2062 

18.8 

109.5 

57.1 

2.5 

2103 

19.0 

110.9 

56.3 

1.7 

2098 

27.0 

77.9 

54.5 

1.7 

2081 

26.8 

77.6 

54.8 

1.7 

2072 

26.7 

77.5 

54.8 

1.1 

2086 

46.9 

44.4 

61.8 

1.2 

2120 

46.0 

46.1 

65.1 

1.2 

2100 

44.3 

47.4 

63.3 

0.6 

2110 

98.5 

21.4 

70.0 

0.6 

2116 

94.4 

22.4 

66.8 

0.6 

2161 

93.3 

23.2 

64.7 

^Compressed  Thickness  Data 
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TABLE  3 
LARGE  SAMPLE  HOLDER  DC  FLOW  RESISTANCE  MEASUREMENT  RESULTS 

(flow  normal  to  surface) 

Compressed  Thickness: 

Mean  R^  =  59.1  x  10 3  MKS  Rayl/m 

Std.  Deviation  =  4.2  x  10 3  MKS  Rayl/m 

Uncompressed  Thickness; 

Mean  R^  =  53.2  x  10  3  MKS  Rayl/m 

Std.  Deviation  =  4.2  x  10 3  MKS  Rayl/m 
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?K  snaller  sannle  holder  was  subsequently  constructed  and  samples  were  cut  which 
allowed  the  measurement  of  flov;  resistance  with  particle  velocity  nomal  to  the 
fiberglass  surface  and  with  particle  velocuty  parallel  to  the  fiberglass 
surface.   Figure  3.3  shows  the  results. 

Mot  only  is  the  DC  flow  resistance  anisotropic  by  a  factor  of  two  to 
three  but  the  effect  of  the  10%  compression  of  the  material  in  the  large  sample 
holder  is  evident.  The  resulting  flow  resistance  parallel  to  the  surface  of 
the  material  was  found  to  be  15,200  ±  500  MKS  Rayl/m  while  that  normal  to  the 
surface  was  34,000  ±  900  PKS  Rayl/m.  Table  4  and  Table  5  give  the  measured 
data. 

2 .  Viscous-Inertial  Transition  Freguency 

The  ratio  of  the  viscous  penetration  depth  to  the  characteristic  pore 
size  is  one  of  the  most  fundamental  dimensionless  parameters  that  determines 
the  nature  of  elastic  wave  propogation  in  a  fluid-saturated  porous  solid.  If 
^visc  <<:  P°re  size,  then  the  motion  of  the  fluid  in  the  pores  is  essentially 
the  same  as  if  the  fluid  were  inviscid.  However,  if  5visc  >>  pore  size  then  the 
fluid  in  the  pores  is  essentially  frozen  to  the  solid,  and  sound  cannot 
propagate  within  the  fluid.  Microscopic  photographs  of  the  fiberglass  material 
provided  an  estimate  of  the  average  distance  between  fibers.  This  mean 
distance  was  determined  to  be  220  ±  80  \m.  The  crossover  from  where  fluid 
inertia  dominates  to  where  flow  resistance  dominates  occurs  when  6visc  ~  1/2 
(mean  pore  diameter).  This  value  was  estimated  to  be  110  ±40  pm.  Since,  the 
thermal  penetration  depth  in  air  is  approximately  the  same  as  the  viscous 
penetration  denth,  the  crossover  from  isothermal  to  adiabatic  compression  in 
the  pores  was  expected  to  be  nearly  coincident  with  the  viscous-inertial 
transition.  Figure  3.4  shows  the  viscous  penetration  depth  over  the 
freguencies  of  interest.    The  transition   freguency  for  this  case  was 


48 


approximate!'/   400    Hz.        Pnr  a  collection  of   straiqht   nines   of      radius   a,    the   IX 
flow  resistance    is  given  bv, 

Rnr  =  ^?  (3.3) 

!JL         a2 

Considering   the  IX!  flow  resitance  measured  with  flow  parallel   to  the 

air-f iberqlass  surface,  the  equivalent  radius  would  he  98   urn.   This  figure 

compares  very  well  with  the  microscopic  observations.  \  radius  of  98  urn  results 
in  a  transition  frequency  of  500  Hz. 

R.   CLASSIC  METHODS  OF  TRANSMISSION  LINE  THEORY 

1 .   Impedance  Tube  Experiment  and  Results 

The  inpedance  tube  technique  is   the  most  widely  used  to  determine  a 

material's  characteristic   impedance,   being  relatively  simple  in  both  theory 

and  measurement.  This  experiment  was  the  first  conducted  of  the  five  classical 

methods. 

\   Rruel  and  Kjaer  Type  4002  Standinq  Have  apparatus  was  used  for  these 

measurements.   The  tube  was  constructed  of  clear  plexiclass  which  allowed 

visual  inspection  of  the  sample.   The  tube  is  1  meter  in  length  and  9.8  cm  in 

diameter.   This  sets  a  lower  limit  on  the  usable  frequency  range  to  be 

approximatelv  85  Hz  for  the  length  of  the  tube  to  be  greater  than  1/4 

wavelength,  and  a  higher  limit  of  approximately  2006  Hz  for  operation  below  the 

lowest  crossmode,  thus  ensuring  only  plane  waves  are  propagated.   the  HP4192A 

Low  Frequency  Impedance  Analyzer  was  used  for  signal  generation  and  detection. 

It  has  a  built  in  freguency  synthesizer  with  a  resolution  of  1  mHz  and  a 

phase-sensitive  detector  with  a  0.1  Hz  bandwidth  to  preclude  noise  and 

harmonics  interfering  with  measurements.   Figure  3.5  shows  the  experimental 

setup.  Since  we  are  primarily  interested  in  the  lowest  longitudnal  mode  in  the 

rectangular  cavity  (particle   velocity  parallel  to  the  surface   of  the 

fiberglass),  which  has  a  freguency  in  the  range  of  600  to  700  Hz,  the  impedance 
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Figure  3.3      Ap     vs     V     in  the  Small  Sample  Holder 
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Figure  3.4  Characteristic  Length  Scale  of  a  Pore  Compared  to  the  Viscous 
Penetration  Depth 
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TABLE   4 


SMALL   SAMPLE   HOLDER  DC    FLOW  RESISTANCE   MEASUREMENT   DATA  AND   RESULTS 

(Flow  Normal   to  Surface) 

Thickness  =  15.60    ±0.01  cm 
Cross-Sectional  Area  =  20.5    ±  0.1   cm2 


Ah 
(cc) 

7.6    ±  0.1 

7.5 

6.0 

5.9 

4.9 

5.0 

4.0 

4.0 

3.0 

3.0 

2.0 

2.0 

1.0 

1.0 


Vol. 
(cc) 

2165    ±  4 

2086 

2052 

2096 

2082 

2108 

2084 

2116 

2150 

2098 

2118 

2151 

2116 

2121 


t 
(sec) 

7.8    ±  0.1 

7.6 

9.0 

9.4 

11.2 

10.9 

13.5 

13.4 

18.3 

18.0 

28.4 

29.1 

57.8 

55.4 


V 

Rdc 

( cc/sec ) 

(MKS 

Rayl/m) 

278 

.6  ±  0.5 

35 

.1> 

:103  ±  9( 

275 

.2 

35 

.1 

229 

.3 

33 

.7 

222 

.3 

34 

.2 

185 

.6 

34 

.0 

192 

.9 

33 

.4 

154 

.3 

33 

.4 

157 

.4 

32 

.7 

117 

.6 

32 

.9 

116 

.6 

33 

.2 

74 

.6 

34 

.6 

74 

.0 

34 

.8 

36 

.6 

35 

.2 

38 

.3 

33 

.7 

Mean     R^  =  34.0   x  10 3  MKS  Rayl/m 
Std  Deviation  =  0.9  x  10 3  MKS  Rayl/m 
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TABLE   5 


SMALL   SAMPLE   HOLDER  DC    FLOW  RESISTANCE   MEASUREMENT   DATA  AND   RESULTS 

(Flow  Parallel   to  Surface) 

Thickness   =   17.80    ±0.01   cm 
Cross-Sectional  Area  =  20.5    ±  0.1   en2 


Ah 

Vol. 

t 

V 

Rdc 

(cc) 

(cc) 

(sec) 

(cc/sec) 

(MKS  Rayl/m) 

3.6  ±  0.1 

2068  ±  4 

7.5  ±  0.1 

276.8  ±  0.5 

14.7x10  3  ±  5( 

3.6 

2100 

7.5 

281.1 

14.4 

3.1 

2066 

9.0 

229.8 

15.2 

3.1 

2120 

9.1 

232.2 

15.1 

2.5 

2109 

11.9 

177.7 

15.9 

2.5 

2124 

11.1 

191.5 

14.7 

1.9 

2126 

14.9 

142.8 

15.0 

1.9 

2115 

15.2 

139.0 

15.4 

1.3 

2075 

22.1 

94.0 

15.6 

1.3 

2108 

22.9 

92.0 

15.9 

0.7 

2150 

42.2 

51.0 

15.5 

0.7 

2104 

39.6 

53.3 

14.9 

Mean     R^  =   15.2   x   10  3  MKS   Rayl/m 
Std.   Deviation  =  0.5   x  10 3  MKS  Rayl/m 
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tube  measurements  were  carried  out  at  frequencies  between  100-1000  Hz  with 
sample  thicknesses  of  4,8,  and  12  inches.  The  different  sample  thicknesses 
were  used  in  order  to  determine  how  thick  the  material  had  to  be  for 
reflections  from  the  riqid  termination  to  be  negligible,  thus  allowing  the 
material's  characteristic  impedance  to  be  determined. 

fleasurements  of  the  reflection  coefficient  from  100  to  1000  Hz  for  all 
three  samples  indicated  that  the  material  could  be  considered  infinite  for 
thicknesses  greater  than  8  to  12  inches.  The  reflection  coefficient,  R,  is 
qiven  by 

R  =  |B/A|  ,  (3.4) 

where  the  ratio  B/A  is  given  by  Eg.  (2.68).  The  reflection  coefficient  changes 
noticeably  in  the  freguency  region  betwen  400  and  500  Hz.  This  is  thought  to 
be  due  to  the  transition  from  isothermal  to  adiabatic  compression  of  the  gas  in 
the  pores. 

The  standing  wave  ratio  and  the  distance  to  the  first  pressure  node 
were  measured  for  the  12"  sample  at  frequencies  between  100  and  1000  Hz.  The 
normalized  characteristic  impedance  lzcnar|/p0c  was  calculated  usinq  the 
"Smith"  chart  nomoqram  with  the  SXJR  and  phase  at  the  first  node  as  input 
parameters.  Table  6  qives  the  results. 

It  was  observed  that  at  frequencies  below  200  Hz  fluctuations  occur  in 
both  maqnitude  and  phase  of  the  impedance.  The  position  and  the  amplitude  of 
the  pressure  antinode  and  node  were  very  difficult  to  measure  accurately  below 
200  Hz,  and  so  these  data  are  not  as  reliable  as  the  data  taken  above  200  Hz. 
From  200  to  400  Hz,  the  impedance  maqnitude  decreases  at  a  larqer  rate  than  at 
frequencies  above  500  Hz,  apparent  evidence  of  the  isothermal  to  adiabatic 
transition  occurrinq  in  the  pores.  This  transition  is  seen  in  both  maqnitude 
and  nhase  as  well  as  the  reflection  coefficient.   The  transition  coincides 
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Figure  3.5     Experimental  Setup  for  the  Impedance  Tube  Method 
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TABLE   6 


IMPEDANCE  TUBE   MEASUREMENT   RESULTS  ON   12"   OF    FIBERGLASS 


Freq. 
(Hz) 

SWR 

Distance 
to  first 
node ( cm ) 

lzcharl/Poc 

^char 
(Deqrees) 

100 

2.79 

59.70 

1.60 

-44.0 

125 

5.83 

59.70 

3.74 

-48.0 

200 

4.79 

36.80 

2.64 

-52.7 

250 

4.40 

28.75 

2.87 

-45.6 

300 

4.00 

24.10 

2.83 

-42.1 

400 

3.46 

18.00 

2.60 

-38.0 

500 

2.23 

14.45 

1.95 

-22.9 

600 

2.47 

11.40 

1.94 

-31.7 

700 

2.26 

9.65 

1.80 

-29.3 

800 

2.30 

8.25 

1.75 

-31.7 

900 

2.15 

7.75 

1.81 

-24.4 

1000 

2.11 

6.90 

1.77 

-25.1 
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with  the  velocity  profile  chanqinq  fron  laminar  flow  in  the  low  frequency  limit 
to  pluq  flow  in  the  hiqh  frequency  limit  (see  Fiqure  2.1). 
2 .   Resonance  Analysis  Experiment  and  Results 

Of  all  the  classic  mehtods,  this  technique  in  most  similar  to  the  new 
technique  in  that  parameters  are  extracted  from  the  resonance  curve  to 
determine  a  sample's  characteristic  impedance.  As  previously  commented,  this 
method  is  limited  to  discrete  frequencies  and  thin  samples  since  the  resonances 
of  the  tube  become  weaker  as  thicker  samples  are  introduced.  The  upper  limit 
thickness  for  our  fiberqlass  was  found  to  be  about  2  inches.  A  1  inch  thick 
sample  was  chosen  since  Yaniv's  method  was  conducted  with  the  same  sample 
thickness  and  since  three  of  the  five  empty  tube  resonances  were  still  distinct 
enouqh  to  make  measurements. 

The  experimental  setup  is  very  similar  to  the  impedance  tube 
experiment.  A  Bruel  and  Kjaer  Type  4002  Standinq  Uave  apparatus  was  used  alonq 
with  a  HP3580A  Spectrum  Analyzer.  The  spectrum  analyzer  was  used  to  determine 
the  resonance  frequencies  and  half-power  widths  of  the  lowest  three  peaks  in 
both  the  empty  tube  and  in  the  tube  filled  with  the  1  inch  thick  fiberqlass 
sample.  A  HP3314A  Function  Generator  was  then  set  to  the  frequency  of  each  of 
these  oeak  and  the  distance  between  the  fiberqlass  surface  and  the  pressure 
maximum  closest  to  the  speaker  was  measured.  This  is  needed  to  calculate  the 
effective  tube  length,  as  an  infinite  impedance  loudspeaker  is  assumed  in  the 
theory.  The  surface  impedance  of  the  sample  was  then  determined  from  these 
measurements  using  Eg.  (2.78).  Figure  3.6  shows  the  experimental  setup  and 
Table  7  gives  the  resulting  values  of  normalized  impedance  values. 

The  lowest  two  resonances  (center  freguencies  of  253  Hz  and  400  Hz, 
respectively)  gave  plausible  values  of  surface  impedance,  however  the  highest 
mode  with  the  absorbing  material  present  (center  frequency  of  538  -  544  Hz)  had 
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Figure  3.6       Equipment  Setup  For  Resonance  Analysis  Method 
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TABLE  7 


RESONANCE  ANALYSIS  SURFACE  IMPEDANCE  RESULTS 
(Impedance  Tube  and  Rectanqular  Cavity) 
(1  Inch  Thick  Fiberglass  Sample) 


Surface  Impedance  Results  in  Impedance  Tube 


Resonance 

Peak 

Number 

Center 
Frequency 
( damped ) 
(Hz) 

HPFW 

(Hz) 

lzsurf l/Poc 

^surf 
( degrees ) 

krr/^n 

1 

253 

23 

2.9 

-85.9 

0.0078 

2 

400 

18 

2.1 

-88.9 

0.0013 

3* 

543 

41 

7.7 

+66.4 

-0.0055 

Surface  Impedance  Results  in  Rectangular  Cavity 


Resonance 
Peak 

Number 


Center 
Frequency 
( damped ) 
(Hz) 


HPFW 

(Hz) 


I z  surf  l/Poc 


i> 


zsurf 
(degrees! 


JSmAia 


841 


126 


1.8 


-65.4 


0.071 


This  point  was  discarded,  as  discussed  in  text, 
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a  higher  Q  than  in  the  empty  tube  and  was  therefore  discarded.  The  erroneous 
data  for  the  third  resonance  peak  is  attributed  to  the  significant  broadening 
of  its  resonance  curve  and  the  resultant  overlapping  with  adjacent  modes. 

This  resonance  analysis  method  was  applied  to  the  rectangular  cavity 
with  and  without  a  1  inch  thick  strip  of  fiberglass  while  exciting  the  [010] 
mode  to  verify  the  results  obtained  in  the  tube.  The  [010]  mode  was  chosen 
because  for  this  mode  the  particle  velocity  is  normal  to  the  fiberglass 
surface,  as  in  the  impedance  tube.  For  these  measurements  in  the  rectangular 
cavity  the  source  was  placed  at  grid  point  (0,0)  (Figure  3.11)  and  the  receiver 
at  grid  point  (20,0).  The  resonance  curve  was  first  measured  in  the  empty 
cavity;  it  had  a  center  freguency  of  873  Hz.  The  normalized  magnitude  of  the 
surface  impedance  obtained  agreed  closely  with  the  value  obtained  by  the 
impedance  tube  for  freguencies  between  800  and  900  Hz.  However  a  significant 
phase  discrepency  was  evident.  The  reason  for  the  discrepency  is  not  known. 
Perhaps  it  is  a  result  of  using  a  point  source  located  in  a  corner  .  Table  7 
lists  the  results. 

As  with  other  impedance  tube  measurements,  this  method  relies  upon  a 
plane  propagating  wave,  where  incidence  on  the  porous  material  is  assumed 
normal  to  the  surface  of  the  sample  under  test.  With  a  point  source  located  in 
one  corner  of  the  cavity  used,  departure  from  plane  wave  propagation  can  be 
expected  with  increasing  sample  thickness  as  the  guality  factor  of  the 
resonance  degrades.  The  extracted  impedances  are  surface  values.  In  both 
geometries  in  which  this  method  was  used  the  magnitude  of  the  surface  impedance 
decreased  with  increasing  freguency  (which  was  expected) .  The  phase  in  all 
measurements  indicates  that  the  impedance  remains  largely  reactive,  and  that 
the  gas  in  the  pores  acts  like  a  stiffness.  Use  of  samples  of  greater 
thickness  would  only  suppress  the  resonance  curves  more. 
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3.   Yaniv's  Method,  Experiment  and  Results 

The  third  and  last  technique  which  uses  transmission  line  theory  is  the 
method  proposed  and  experimentally  verified  by  Yaniv  [Ref .  5] .  The  experiment 
was  conducted  using  the  impedance  tube  of  the  previous  two  experiments,  with 
several  modifications.  Yaniv's  method  requires  measurinq  the  surface  impedance 
with  the  sample  backed  by  a  riqid  termination  and  by  an  acoustic  short-circuit, 
obtained  by  movinq  the  riqid  termination  such  that  a  1/4  X  air  gap  exists 
between  the  rear  surface  of  the  material  and  the  riqid  termination. 

The  riqid  end  cap  used  in  the  two  previous  experiments  was  replaced 
with  a  2  inch  thick  disc  of  plexiqlass  with  a  3/8  inch  diameter  hollow  aluminum 
tube  threaded  into  the  side  away  from  the  driver.  A  small  amount  of  hiqh 
vacuum  qrease  was  used  to  seal  the  edqes  of  the  disc  at  the  tube  walls  while 
still  allowinq  for  positioning.  The  front  face  of  the  disc  had  a  recess  cut 
out  of  it  into  which  a  FET  microphone  mounted  at  the  end  of  a  3/16  inch  hollow 
aluminum  rod  could  fit  snuqly  and  flush  so  as  to  provide  a  smooth  surface  when 
the  disc  was  placed  aqainst  the  back  end  of  the  sample.  The  3/16  inch  diameter 
rod  fit  inside  the  3/3  inch  diameter  rod,  providinq  support  and  allowinq  the 
microphone  to  be  positioned  anywhere  between  the  backside  of  the  porous  solid 
and  the  riqid  termination.  Electrical  connection  to  a  power  supply  switchbox 
was  also  made  by  wires  which  passed  throuqh  the  small  hollow  rod.  This 
conf iquration  allowed  the  disc  and  microphone  to  serve  as  the  riqid  termination 
for  the  open  circuit  impedance  measurement  and  also  allowed  accurate 
positioninq  of  the  termination  to  provide  the  required  one-quarter  wavelength 
air  gap  between  the  porous  solid  and  the  rigid  termination  at  each  freguency. 
Figure  3.7  shows  the  movable  termination  and  Figure  3.0  shows  the  eguipment  set 
setup. 

A.  1  inch  thick  disc  of  fiberglass  was  inserted  into  the  impedance  tube 
such  that  the  distance  betwen  the  driver  and  the  sample  was  at  least  a  guarter 
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Figure  3.7  The  Rigid  Sliding  Termination  (Yaniv'S  Method) 


wavelength.  The  lowest  usable  freguency  was  therefore  approximately  250  Hz, 
the  upper  limit  remaining  approximately  2000  Hz  to  preclude  the  lowest 
crossmode.  Measurements  were  made  of  the  first  minimum  and  maximum  pressure 
amplitudes  in  front  of  the  sample  and  of  the  distance  between  the  sample 
surface  and  the  first  pressure  minimum  in  from  of  the  sample.  For  each 
freguency,  these  measurements  were  performed  under  conditions  of  open  and  short 
circuit  terminations.  A  short  circuit  termination  at  the  back  side  of  the 
fiberglass  disc  was  obtained  by  adjusting  the  position  of  the  plexiglass  plug 
until  a  null  output  was  observed  from  the  microphone,  which  was  held  flush  to 
the  fiberglass  disc.  It  was  found  that,  within  experimental  error,  the  result- 
ing distance  between  the  porous  disc  and  the  plexiglass  plug  was  within  1% 
of  the  1/4  wavelength  determination  made  in  the  empty  tube  at  the  freguencies 
used  in  the  experiment.  The  standing  wave  ratio  and  the  reflection  coefficient 

phase  angle  were  calculated  using  Eg.  (2.67)  and  Eg.  (2.70),  respectively.  The 
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Figure  3.8  Equipment  Setup  for  Yaniv's  Method 
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open  and  short  circuit  impedances  were  found  using  the  "Smith"  chart.  The 
characteristic  impedance  was  then  found  by  use  of  Eg.  (2.84). 

Determination  of  the  complex  propagation  constant,  y  =  jk,  reguired 
finding  the  roots  of  Eg.  (2.85)  such  that, 

tanh[yL]  -  [Zss/?os]  l  /2   =  0   .  (3.5) 

The  roots  were  found  using  a  numerical  technigue  known  as  Ward's  method.  The 
details  are  given  in  Appendix  A  [Ref .  19] .  The  results  of  this  experiment  are 
listed  in  Table  8. 

As  for  the  previous  methods,  a  decrease  in  impedance  magnitude  with 
increasing  freguency  was  observed,  with  the  transition  from  isothermal  to 
adiabatic  compression  between  400  to  700  Hz  in  this  1  inch  thick  sample.  The 
complex  wave  number  also  shows  this  transition.  The  data  and  results  for  400 
Hz  appear  to  anomalous.  The  reason  is  not  known. 

C.   PROBE  TUBE  METHOD,  EXPERIMENT  AND  RESULTS 

The  probe  tube  measuremnets  were  made  by  setting  up  a  plane  wave  in  the 
Bruel  and  Kjaer  Type  4002  Standing  Wave  apparatus  filled  partially   with 

fiberglass  and  exploring  the  sound  dield  by  means  of  a  probe  tube  extension  to 

the  impedance  tube  microphone.   The  slope  of  the  decay  of  the  pressure 

amplitude  with  distance  provided  a  measurement  of  the  attenuation  constant  a. 

The  change  in  phase  of  the  oscillatory  wave  with  distance  provided  the  wave 

/v 

number  k.  The  complex  wave  number  is  then  found  from  k  =  k  -  ja. 

The  equipment  setup  is  shown  in  Figure  3.9.  Three  specimens,  of 
thickness  4,  8,  and  12  inches,  were  chosen  for  measurement  and  comparison  with 
the  impedance  tube  results.  Each  sample  thickness  was  made  up  of  the  reguired 
number  of  fiberglass  discs,  9.84  cm  diameter  and  4.76  cm  thick,  inserted  into 
the  tube  opposite  the  driver,  taking  care  to  avoid  compression  of  each  disc 
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TABLE   8 


YANIV'S   METHOD   EXPERIMENTAL   RESULTS 


JOS 


'SS 


Freq. 
(Hz) 

S^JR 

Distance  to 
first  node 
(cm) 

s\m 

Distance  to 
first  node 
(cm) 

250 

50.40 

30.90 

1.36 

30, 

.50 

300 

42.40 

25.80 

1.33 

25, 

.70 

400 

24.03 

18.60 

1.89 

33, 

.20 

500 

24.71 

14.10 

1.68 

21, 

.70 

600 

12.29 

11.50 

1.65 

17, 

.20 

700 

9.57 

9.70 

1.56 

14 

.70 

800 

8.71 

8.20 

1.70 

12 

.60 

900 

7.16 

6.80 

1.70 

11 

.15 

1000 

6.42 

5.90 

1.67 

10 

.45 

Freq. 
(Hz) 

lzcharl/P0c 

^char 
( degrees ) 

k 
(Rad/m) 

a 

(Np/m) 

cph 
(m/s) 

250 

2 

.79 

-44.8 

12.82 

13, 

.96 

122.5 

300 

2 

.78 

-43.6 

12.56 

13, 

.86 

150.2 

400 

2 

.04 

-22.3 

15.17 

8, 

.80 

165.7 

500 

2 

.17 

-30.1 

20.82 

12, 

.75 

150.9 

600 

2 

.10 

-28.9 

21.70 

15 

.77 

173.7 

700 

1 

.98 

-28.4 

21.73 

16 

.81 

202.2 

800 

1 

.94 

-27.9 

24.88 

17 

.50 

202.1 

900 

1 

.75 

-27.7 

27.75 

18 

.40 

203.7 

1000 

1 

.62 

-25.6 

28.41 

17 

.92 

221.2 
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while  trying  to  minimize  any  air  gaps  between  disc  surfaces.  Each  disc  had  a 
0.71  cm  diameter  hole  cut  in  the  center  for  a  close  fit  around  the  probe  tube 
without  any  deformation  of  the  fiberglass. 

The  probe  tube  was  manufactured  out  of  steel  tubing  with  the  same  outside 
and  inside  diameter  of  the  existing  open  ended  probe  tube  on  the  standing  wave 
apparatus.  A  copper  cap  was  welded  onto  the  end  of  the  extension  next  to  a 
3/16  inch  diameter  hole  to  allow  the  measurement  of  the  sound  field.  The 
extension  attached  to  the  existing  probe  tube  with  a  plastic  sleeve.  This  is 
shown  in  Figure  3.10. 

The  experiment  was  conducted  at  freguencies  of  125,  250,  500,  and  1000  Hz 
for  each  of  the  three  sample  thicknesses,  but  only  the  12  inch  thick  sample 
was  long  enough  so  that  the  wave  reflected  from  the  rigid  end  of  the  tube  was 
sufficiently  attenuated  to  be  negligible  compared  to  the  incident  wave.  The 
pressure  amplitude  at  the  rigid  termination  was  measured  to  be  20  dB  less  than 
that  at  the  sample  boundary.  Determination  of  the  attenuation  constant  was 
made  by  plotting  pressure  versus  distance  on  semi-log  paper  and  extracting  the 
slope.  The  phase  of  the  pressure  was  measured  with  a  SR-350  Dual  Channel 
Lock-In  Amplifier.  The  wave  number  was  determined  from  the  change  of  phase 
with  distance  into  the  fiberglass.  Table  9  displays  the  experimental  results 
in  terms  of  k  and  a. 

The  impedance  magnitude  decreases  with  freguency  while  the  phase  indicates 
that  the  material  becomes  more  resistive  with  increasing  frequency.  Only  the 
absorption  coefficient  shows  any  sign  of  the  transition  between  isothermal  and 
adiabatic  compressions  in  the  pore  volume. 

D.   MATERIAL  PROPERTY  METHOD,  EXPERIMENT  AND  RESULTS 

All  the  methods  discussed  so  far  involve  measuring  certain  characteristics 
of  a  progressive  plane  wave  in  the  porous  material  or  its  proximity.   Beranek 
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Figure  3.9  Equipment  Setup  for  Probe  Tube  Experiment 
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has  shown  that  the  complex  wave  number  and  characteristic  impedance  can  also  be 
derived  from  the  basic  physical  properties  of  the  material  with  the  aid  of 
charts.  These  properties  are: 

a)  specific  flow  resistance  per  unit  length,  R^(~  R^) 

b)  porosity,  Y 

c)  structure  factor,  kst  (tortuosity) 

d)  volume  coefficient  of  elasticity  of  air  in  the  interstices  K  (Bulk 
modulus) 

e)  volume  coefficient  of  elasticity  of  the  skeleton  of  the  material. 
Measurements  of  the  specific  flow  resistance  and  porosity  have  already  been 
discussed.   The  remaining  three  material  properties  are  difficult  if  not 
impossible  to  measure,  and  will,  based  on  a  reasonable  knowledge  of  the 
material,  be  approximated. 

The  structure  factor  kst,  known  as  the  tortuosity  in  Biot  theory,  is 
essentially  a  path  length  correction  factor.  Experimentally,  it  is  found  to 
depend  most  strongly  on  the  porosity  of  the  porous  solid  [Ref .  11] .  In  the 
limit  porosity  =  1,  the  tortuosity  also  eguals  1  for  typical  porous  maerials, 
and  certainly  for  fiberglass.  As  the  porosity  of  the  samples  used  in  this 
research  was  found  to  be  very  close  to  1,  a  value  of  1  for  the  tortuosity  was 
assumed  in  all  calculations. 

The  volume  coefficient  of  elasticity  of  air  in  the  interstices,  K,  is  an 
important  term;  if  isothermal  conditions  can  be  assumed  then  K  =  1  x  10 6  dynes/ 
cm2;  if  adiabatic  conditions  exist  then  K  5  1.4  x  106  dynes/cm2.  For  most 
fibrous  acoustic  materials  however,  K  will,  in  seme  freguency  region,  be 
neither  isothermal  nor  adiabatic  but  will  be  something  in  between  and  will  be 
complex  valued  and  freguency  dependent,  as  is  our  case. 

Finally,  the  volume  coefficient  of  elasticity  of  the  skeleton,  K^g^/ 
although  not  directly  involved  in  the  calculation  of  k  and  Zchar'  nust  be 
known  or  approximated  since  it  is  the  ratio  of  K/Kg^gi  that  determined  whether 

67 


c 
o 

•H 
CO 

c 
o 

— I 

X 

w 

I 


CD 
Li 

01 


68 


TABLE   9 


PROBE  TUBE   EXPERIMENTAL   RESULTS 
12"   Thick   Fiberqlass  Sample) 


Freq 
(Hz) 

k 
(Rad/n) 

a 

(Np/m) 

125 

14.8 

13.6 

250 

20.9 

16.3 

500 

29.7 

22.7 

1000 

43.3 

21.7 

cph 
(n/s) 

53.1 

75.2 

105.8 

145.1 
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a  material  is  treated  as  a  soft  blanket  or  a  rigid/dence  blanket.  Based  on  the 
appearance,  texture,  and  high  flow  resistivity  of  the  fiberglass  it  was  assumed 
that  the  material  should  be  treated  as  a  rigid/dense  blanket.  [Ref.  4,  pp. 
839-840] 

Once  the  properties  of  the  material  were  characterized,  Eg.  (2.89)  could  be 
used  to  directly  calculate  or  Figure  2.6  and  Figure  2.7  used  to  graphically 
determine  the  magnitude  and  phase  of  the  complex  propagation  constant  b  e  jk  as 
a  function  of  freguency,  from  which  the  complex  wave  number  is  determined. 
Once  these  were  determined  the  characteristic  impedance  could  be  calculated  by 
use  of  Eg.  (2.90).  The  values  determined  from  the  graphs  for  k  and  the 
calculated  values  of  Zchar  assuming  isothermal  conditions  in  the  pore  volume 
are  listed  in  Table  10.  The  same  guantities,  assuming  adiabatic  conditions, 
are  listed  in  Table  11. 

As  expected,  the  impedance  magnitude  decreases  with  increasing  freguency. 
The  impedance  phase  calculated  for  both  isothermal  and  adiabatic  conditions  is 
the  same,  a  conseguence  of  the  bulk  modulus  being  a  real  guantity  in  the  two 
limiting  cases.  If  isothermal  conditions  were  to  exist,  then  the  impedance 
phase  would  approach  an  angle  of  -45  degrees  (pure  diffusion  of  sound 
pressure).  Since  the  calculated  phase  values  remain  nearer  to  the  diffusive 
limit,  the  assumption  that  isothermal  conditions  exist  in  the  pore  volume  would 
be  most  nearly  correct. 

Over  the  frequency  range  considered,  the  isothermal  values  for  k  and  a  are 
greater  than  the  adiabatic  values,  which  is  consistent  with  the  fact  that  the 
isothermal  sound  speed  is  less  than  the  adiabatic  sound  speed. 

E.   THE  NEW  TECHNIQUE,  EXPERIMENT  AND  RESULTS 

The  plane  wave  resonator  chosen  to  conduct  all  experiments  was  constructed 
out  of  plexiglass  with  interior  dimensions  as  follows: 
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Lenqth  (x-direction)  =  25.4  ±0.1  cm 

Width  (y-direction)  =  19.3  ±  0.1  en 

Heiqht  ( z-direction)  =  7.3  ±  0.1  cm. 
This  particular  cavity  was  constructed  for  a  laboratory  experinent  on  standing 
waves  in  a  rectangular  room.  The  minimum  wall  thickness  was  1.27  cm  (1/2"). 
The  cavity  was  mounted  on  three  legs  with  the  source  hole  located  in  the  floor 
at  one  corner  to  allow  the  excitation  of  all  modes.  The  top  cover  of  the  cavity 
is  removable  to  allow  insertion  of  fiberglass  samples.  When  installed,  the  ton 
cover  seals  the  cavity  and  slides  to  allow  the  receiver  position  to  be  varied. 
A  small  microphone  receptacle  mounted  on  a  slidebar  installed  in  the 
removable  cover  allows  the  pressure  to  be  measured  throughout  the  cavity  with  a 
probe  tube.  All  sliding  surfaces  were  constructed  with  very  close  tolerances 
and  coated  with  high  vacuum  grease  to  reduce  leakage  and  to  provide  a  tight 
fit  when  the  cavity  cover  was  installed. 

The  receiver  slide  bar  assembly  consisted  of  three  portions:  the  slide  bar, 
the  microphone  receptacle,  and  the  receiver  microphone.  The  slide  bar  fit  into 
a  groove  machined  into  the  interior  surface  of  the  top  cover  and  fit  flush  with 
this  surface.  A  threaded  hole  was  cut  into  the  slidebar  to  allow  mounting  the 
microphone  receptacle.  The  microphone  receptacle  consisted  of  a  hollow 
aluminum  cylinder  with  a  thumb  screw  to  set  and  hold  the  height  of  a  brass 
probe  tube  assembly.  The  probe  tube  assembly  allowed  insertion  of  a  General 
Radio  precision  microphone  flush  with  the  smaller  bore  probe  tube  extension 
which  could  be  inserted  into  the  cavity  interior  without  disrupting  the  sound 
field. 

The  bottom  of  the  cavity  was  marked  with  grid  lines,  21  in  the  length 
direction,  15  in  the  width  direction,  spaced  1/2"  apart,  producing  336 
individual  grid  points.   The  in-phase  and  guadrature  components  of  pressure 


71 


TABLE    10 


MATERIAL   PROPERTIES   METHOD   RESULTS 

( I so the  rma 1  Cond  i  t  ions ) 

(RqC  =   15,000   MFCS  Rayl/m) 


Freq. 
(Hz) 

|2Cnarl/P0c 

^char 
( degrees ) 

k 

(Rad/m) 

a 
(Np/m) 

cph 
(m/s) 

100 

3.66 

-43.5 

6.83 

6.48 

92.0 

125 

3.17 

-43.0 

7.47 

6.96 

105.1 

200 

2.64 

-42.0 

10.08 

9.08 

124.7 

250 

2.32 

-41.5 

11.17 

9.89 

140.6 

300 

2.12 

-41.0 

12.33 

10.72 

152.9 

400 

1.83 

-39.0 

14.66 

11.82 

171.4 

500 

1.66 

-37.5 

16.94 

13.00 

185.5 

600 

1.49 

-37.0 

18.37 

13.84 

205.2 

700 

1.40 

-35.0 

21.26 

14.88 

206.9 

800 

1.32 

-33.5 

22.63 

14.98 

222.1 

900 

1.24 

-32.0 

24.46 

15.28 

231.2 

1000 

1.20 

-31.0 

26.38 

15.85 

238.2 
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TABLE  11 


MATERIAL  PROPERTIES  METHOD  RESULTS 
(Adiabatic  Conditions) 
(R^  =  15,000  MKS  Rayl/m) 


Freq. 
(Hz) 

lzchar!/P0c 

zchar 
( degrees ) 

k 
(Rad/m) 

(Np/m) 

°ph 
(m/s) 

100 

4.28 

-43.5 

5.69 

5.40 

110.4 

125 

3.93 

-43.0 

6.60 

6.16 

119.0 

200 

3.04 

-42.0 

8.31 

7.48 

151.2 

250 

2.70 

-41.5 

9.29 

8.22 

169.1 

300 

2.39 

-41.0 

9.96 

8.65 

189.3 

400 

2.09 

-39.0 

11.91 

9.65 

211.0 

500 

1.88 

-37.5 

13.71 

10.52 

229.1 

600 

1.71 

-37.0 

15.05 

11.34 

250.5 

700 

1.68 

-35.0 

17.65 

12.36 

249.2 

800 

1.54 

-33.5 

18.86 

12.48 

266.5 

900 

1.44 

-32.0 

20.14 

12.59 

280.8 

1000 

1.39 

-31.0 

21.82 

13.11 

287.9 
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were  measured  at  each  grid  point  for  the  node  and  sample  selected.  Fiqure  3.11 
shows  the  location  of  the  source  and  the  grid  numbering  schene  that  was  applied 
to  the  interior  floor  of  the  cavity. 

Ten  fiberglass  samples  were  cut  from  a  fiberglass  sheet,  with  thicknesses 
ranqing  fron  0.48  cm  (3/16  inch)  to  4.76  cm  (1  7/8  inch)  the  maximum  thickness 
being  egual  to  the  full  thickness  of  the  sheet.  The  samples  were  cut  in  the 
form  of  rectangular  strips  so  as  to  fit  snugly  along  the  side  wall  opposite  the 
source.  Figure  3.11  shows  the  location.  The  strips  extended  from  grid  point 
(0,15)  to  grid  point  (20,15).  The  samoles  were  cut  such  that  the  surface 
normal  to  the  flow  velocity  in  the  original  flow  resistance  experiment  would 
also  be  normal  to  the  particle  velocity  expected  when  the  [010]  mode  was 
excited  (i.e.,  the  "width"  direction  in  Figure  3.11)  and  parallel  to  the 
particle  velocity  when  the  [100]  mode  was  excited  (the  "length"  direction). 
Table  12  lists  the  samples  used  and  their  thicknesses. 

Three  experiments  were  conducted: 

1 .  The  mode  structure  in  the  cavity  was  measured  over  the  freguency  range 
of  interest  for  each  sample.  The  effect  of  introducing  the  fiberglass 
samples  on  the  strength  and  shape  of  the  resonances  was  noted.  Of 
particular  interest  were  the  [100]  and  [010]  modes,  which  were  expected 
to  be  relatively  distinct. 

2.  The  DC  flow  resistance  of  the  air-filled  fiberglass  was  then  extracted 
from  the  above  measurements  using  the  dispersion  relation  developed 
in  Chapter  2  for  the  allowed  modes  (Eg.  (2.51)),  following  which  the 
complex  wave  number  and  characteristic  impedance  were  calculated. 

3.  The  pressure  field  of  the  [100]  and  [010]  modes  was  measured  for 
varying  sample  thicknesses  and  compared  to  that  predicted  by  theory. 

The  equipment  setup  used  to  conduct  all  three  of  these  experiments  is  shown 
in  Figure  3.12. 
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Figure  3.11  Rectangular  Cavity  Grid  Map 

1 .  Modal  Freguencies  Up  to  2000  Hz  for  Various  Fiberglass  Thicknesses 

The  first  of  the  three  experiments  was  conducted  to  determine  the  modes 
present  in  the  empty  cavity  and  the  effects  of  increasing  sample  thickness  on 
the  shape  and  strength  of  each  resonant  mode.  From  mode  theory,  seven  modes 
were  expected  to  exist  in  the  empty  cavity  for  frequencies  less  than  2000  Hz. 
Of  these  seven  modes,  the  two  lowest  were  the  [100]  and  [010]  modes,  both  of 
which  occur  below  1000  Hz.  All  seven  modes  were  observed,  with  the  two  lowest 
modes  being  strong  and  relatively  distinct  compared  to  the  others.  This  is 
shown  in  Figure  3.13.  Table  13  lists  the  theoretically  calculated  and 
experimentally  measured  peak  freguencies  of  the  modes  present  in  the  cavity  up 
to  2000  Hz. 
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There  was  some  concern  that  as  thicker  samples  were  introduced  the 
resonance  peaks  would  overlap,  makinq  the  measurement  of  the  half-power  widths 
unreliable.  This  condition  set  the  upper  limit  on  usable  sample  thickness  to  1 
7/8"  for  the  [100]  mode  (sample  #10)  and  11/16"  for  the  [010]  mode  (sample 
#5). 

For  all  observed  modes,  the  center  frequency  down-shifted  as  the  sample 
thickness  was  increased,  while  the  0  decreased.  The  [100]  mode  has  an  empty 
cavity  peak  frequency  and  Q  of  680.2  Hz  and  147.9,  respectively.  Insertinq  the 
1  7/8"  (sample  #10)  fiberqlass  strip  down-shifted  its  center  frequency  and 
lowered  the  Q  to  570.0  Hz  and  3.2,  respectively,  for  a  total  frequency  chanqe 
of  110.2  Hz.  The  [010]  mode  has  a  center  frequency  of  870.2  Hz  and  a  0  of  57.6 
in  the  empty  cavity,  but  when  the  11/16"  (sample  #5)  fiberqlass  strip  was 
inserted  the  peak  frequency  only  down-shifted  17  Hz  to  853.2  Hz  with  a  Q  of 
21.1.  For  both  mode  measurements  the  receiver  probe  tube  remained  flush  with 
the  ceilinq  of  the  cavity  at  qrid  point  (20,0).  See  Fiqures  3.14  and  3.15 

Several  observations  from  this  experiment  were  made: 

(1)  The  Q  and  the  center  frequency  of  both  modes  decreased  with  increasinq 
sample  thickness. 

(2)  Compared  to  the  [010]  mode,  the  [100]  mode,  with  particle  velocity 
parallel  to  the  fiberqlass  surface,  allowed  thicker  samples  to  be 
introduced  before  the  resonance  curve  spread  such  that  the  Q  could  no 
longer  be  reliably  determined. 

(3)  For  all  sample  thicknesses,  the  Q  of  the  resonance  curves  was 

larger  with  the  receiver  positioned  above  the  source  (grid  point  (0,0)) 

than  in  the  opposite  corner  (grid  point  (20,0)).   This  was  assumed 

to  be  due  to  near  field  effects  of  the  source  on  the  receiver. 

Consequently,  only  the  resonance  data  taken  with  the  receiver  at  qrid 

point  (20,0)  is  reported  here. 
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TABLE  12 


FIBERGLASS  SAMPLE  THICKNESS 


Thickness 

Thickness 

Sanple  # 

( inches) 

(cm) 

1 

3/16  ±  1/32 

0.48  ±  0.01 

2 

1/4 

0.64 

3 

3/8 

0.95 

4 

1/2 

1.27 

5 

11/16 

1.75 

6 

7/8 

2.22 

7 

1 

2.54 

8 

1  1/4 

3.18 

9 

1  1/2 

3.81 

10 

1  7/8 

4.76 
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Figure  3.12  Experimental  Setup  for  the  New  Method 
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TABLE  13 


RESONANCE  MODES  OF  THE  EMPTY  CAVITY 
(Theoretical  and  Measured) 


Mode  f 

100 

010 

110 

200 

210 

020 

120 


Resonant 

Frequency 

(calculated) 

679  ±  4 

871  i  5 

1105  ±  9 

1358  ±  12 

1614  1   14 

1742  ±  19 

1870  ±  23 


Resonant 

Frequency 

(measured) 

676  ±  1 

868 
1099 
1354 
1608 
1738 
1865 
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Fiqure  3.13  Mode  Structure  of  the  Empty  Cavity 
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(4)   Overlap  of   the   [110]   mode  with  the   [010]   mode  was  evident  for 

fiberglass  thickness  of  1/2"  (sample  #4)   and  greater.   Overlap  of 

the  [100]   mode  with  the   [010]   mode  was  not  significant  until  a 

thickness  of  1  1/4"  (sample  #8)  of  fiberglass  was  introduced.  See 

Figures  3.14  and  3.15. 

2.  Measured  fc  and  Q  of  the  [100]  and  [010]  Modes  for  Various  Fiberglass 
Thicknesses 

The  second  experiment  consisted  of  measuring  the  center  freguency  and  Q 
of  the  [100]  and  [010]  modes  as  a  function  of  sample  thickness.  The  receiver 
was  placed  in  the  ceiling  of  the  cavity  at  grid  point  (20,0)  at  the  opposite 
end  of  the  long  wall  from  the  source.  This  location  was  chosen  to  minimize 
direct  acoustic  cross-talk  between  the  source  and  receiver.  For  the  [100] 
mode,  each  of  the  10  samples  was  inserted  and  the  center  freguency  and  0  were 
measured  from  the  resonance  curve.  This  sample  procedure  was  repeated 
for  the  [010]  mode  except  only  samples  1  through  5  were  used.  The  center 
freguency  and  half-power  freguencies  were  mesured  to  within  a  tenth  of  a  Hertz 
by  reading  the  freguency  of  the  1  to  1.5  MHz  Local  Oscillator  output  of  the 
HP3580A  Spectrum  Analyzer  on  a  HP5300B  Freguency  Counter.  Setting  the  scale  to 
1  dB/div  on  the  Spectrum  Analyzer  allowed  the  center  freguency  and  half-power 
full  width  of  the  resonance  curves  to  be  measured.  Tables  14  and  15  list  the 
experimentally  determined  values  of  fc  and  Q  for  the  [100]  and  [010]  modes. 

Following  the  experimental  measurements,  theoretical  values  of  fc  and 
0  for  the  [100]  and  [010]  modes  as  a  function  of  sample  thickness  were 
calculated  from  the  complex  roots  of  the  dispersion  relation  (Eg.  (2.51))  for 
various  values  of  DC  flow  resistance  using  a  numerical  analysis  technigue  known 
as  Ward's  method  [Ref.  19,  pp.  169-170].  A  brief  description  of  this  method  is 
given  in  Appendix  A.   However,  before  a  comparison  between  theory  and 
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Figure  3.14     Mode  Structure  of  Cavity  With  a  11/16"   Fiberglass  Strip 
(sample  #5) 
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Figure  3.15  Mode  Structure  of  the  Cavity  With  1  7/8"  Fiberglass  Strip 
(sample  #10) 
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experiment  could  be  made,  and  the  best  fitting  value  of  DC  flow  resistance 
determined,  it  was  necessary  to  correct  one  or  the  other  for  empty  cavity 
losses.  Possible  sources  of  empty  cavity  losses  are  considered  in  the  next 
subsection. 

3.  Sources  of  Loss  in  the  Empty  Cavity 

In  derivinq  the  dispersion  relation  for  the  new  resonance  technique  it 
was  assumed  that  the  only  losses  which  occur  in  the  plane  wave  resonator  are 
due  to  viscous  and  thermal  losses  within  the  fluid  saturated  porous  solid. 
Althouqh  these  are  the  larqest  sources  of  loss  in  the  cavity,  other  los 
mechanisms  must  be  accounted  for.  To  account  for  these,  several  assumptions 
were  made: 

(1)  the  losses  in  the  empty  cavity  remain  the  same  even  thouqh  a 
portion  of  the  volume  and  wall  surface  is  lost  due  to  the  introduction 
of  a  fiberqlass  sample  (less  than  10%), 

(2)  the  shift  of  resonance  frequency  due  to  the  introduction  of  a  fiber- 
qlass sample  does  not  effect  the  Q  of  the  empty  cavity  and, 

(3)  each  loss  mechanism  can  be  treated  as  actinq  independently. 

Generally,  stated,  it  was  assumed  that, 

Total  losses  in  the  =  losses  in  the  +  losses  due  to 
fiberqlass   filled  empty  cavity  only  to  fiber- 

cavity  qlass 

which  can  be  written  in  terms  of  quality  factors  as 

1 „!_+!_  (3#6) 

Otgt  Oe  °f 

where  Q^cr  is  the  0  of  the  cavity  due  to  all  loss  mechanisms,  Qj?  is  the  empty 
cavity  Q  and  Op  is  the  0  due  to  viscous  and  thermal  losses  within  the 
fiberqlass.  Eq.  (3.6)  was  used  to  compare  the  measured  values  of  Opcr  to  the 
calculated  values  of  Op  for  the  [100]  and  [010]  modes  by  settinq  the  values  of 
Qp  equal  to  those  experimentally  determined  for  the  empty  cavity. 
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TABLE    14 


[100]    PDDE   CENTER   FREQUENCY  AND     Q     VALUES 


Thickness 

^c 

E* 

fu 

Sanple  # 

(cm) 

(Hz) 

(Hz) 

(Hz) 

Q 

Empty 

0.0 

680.2  ±  0.1 

678.2  ±  0.1 

682.8  ±  0.1 

147.9  i  0. 

1 

0.48 

669.5 

663.9 

675.3 

58.7 

2 

0.64 

667.5 

661.3 

674.1 

52.2 

3 

0.95 

658.6 

648.8 

668.1 

34.1 

4 

1.27 

652.1 

640.1 

664.2 

27.1 

5 

1.75 

646.3 

628.9 

659.9 

20.9 

6 

2.22 

623.6 

587.0 

650.5 

9.8 

7 

2.54 

614.1 

567.4 

647.6 

7.7 

8 

3.18 

592.3 

531.9 

629.7 

6.1 

9 

3.81 

571.6 

506.6 

624.3 

4.9 

10 

4.76 

570.0 

495.4 

672.6 

3.2 
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TABLE    15 


[010]    fDDE   CENTER   FREQUENCY   AND     0     VALUES 


Sample  # 

Thick 
(cm) 

ness 

Ec 
(Hz) 

(Hz) 

(Hz) 

0 

Empty 

0.0 

870.2 

±   0.1 

863.0 

±  0.1 

878.1  t 

0.1 

57.6 

t  0.9 

1 

0.48 

865.1 

854.5 

874.4 

43.5 

2 

0.64 

864.4 

852.0 

876.2 

35.7 

3 

0.95 

858.8 

846.4 

872.3 

33.2 

4 

1.27 

856.5 

841.1 

873.1 

26.8 

5 

1.75 

853.2 

834.4 

874.9 

21.1 
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An  attempt  to  identify  the  major  source  of  loss  in  the  empty  cavity  was 
made,  as  this  in  turn  determined  the  value  of  Qg  to  be  used  in  correct inq  for 
losses  other  than  in  the  fiberglass.  The  loss  mechanisms  in  the  empty  cavity 
were  assumed  to  be  due  to  3  mechanisms: 

—  transmission  of  enerqy  through  the  plexiglass  walls  of  the  cavity, 

—  viscous  and  thermal  losses  at  the  interior  surfaces  of  the  cavity, 

—  leakage  of  energy  due  to  incomplete  sealing  of  the  top  cover  and  sliding 
section  of  the  cavity  and  slide  bar  assembly. 

The  transmission  of  energy  through  the  walls  was  of  concern  mainly  since  the 

cavity  was  not  designed  specifically  for   this  experiment.   For   normal 

incidence,  the  transmission  loss  through  the  walls  of  the  cavity  was  calculated 

to  be  32  dB  for  the  minimum  wall  thickness  (1.27  cm),  [Ref.  17,  p.  128]  .  Based 

on  the  maximum  wall  thickness  (2.54  cm)  the  lowest  expected  coincidence 

freguency  for  a  free  plexiglass  plate  was  calculated  to  be  1282  Hz.   For  a 

plate  fixed  about  its  edges  this  frequency  is  higher.   Therefore,  transmission 

by  coupling  to  flexural  waves  in  the  walls  can  be  ignored.   The  following 

eguation  was  used  to  estimate  the  0  resulting  from  each  source  of  loss: 

2  it  Estorecj 

0  =  f; — ; — T"  '  (3*7) 

Llost/cycle 
where 

^stored  =  eV  (3'8) 

and  e  is  the  average  energy  density,  and  V  is  the  interior  volume  of  the 

cavity.  The  energy  lost  per  cycle  can  be  written 

Elost/cycle  "  1*  e  TI  f  '  (3,9) 

where  A  is  the  surface  area,  c  is  the  speed  of  sound  in  the  fluid,  Tj  the 

power  transmission  coefficient  and  f  is  the  frequency.   Inserting  Eq.   (3.8) 

and  (3.9)  into  Eq.  (3.7)  results  in 
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O  —X7 

Otrans  =  "^  •  (3.10) 

For  loss  due  to  transmission  throuqht  the  cavity  walls,  assuming  1000 

Hz  as  the  worst  case  (A  =  0.343  m) ,  the  value  of  Qtrans  was  calculated  to  be 

2700.    It  was  concluded  that  transmission  through  the   walls  was  not 

significant. 

Next  considered  were  the  effects  of  the  viscous  and  thermal  losses 
occurring  at  the  cavity  surfaces.  Again  Eg.  (3.7)  was  used  to  determine  QSurf' 
with  Est-ore(j  remaining  the  same  as  in  Eg.  (3.8).  The  energy  lost  per 
cycle  due  to  viscous  effects  can  be  written  approximately  as 

Elost/cycle  -     ^  5visc  (3.11) 

where  <$visc  *s  t^ie  viscous  penetration  depth  in  air.  The  energy  lost  per 
cycle  due  to  thermal  effect  can  be  written  approximately  as 

E lost/cycle  2  eA  (rD«th  (3.12) 

where  y  is  the  ratio  of  specific  heats  Cp/cv  (=  1.4)  for  the  air  and  5^  is 
the  thermal  penetration  depths  in  air.  Since  the  viscous  and  thermal 
penetration  depths  in  air  are  approximately  egual,  an  expression  for  the 
energy  lost  per  cycle  due  to  both  viscous  and  thermal  effects  is 

Elost/cycle  -     A   'Vise  Y  •  (3.13) 

In  other  words, 

Svisc  +  (^1)(Sth  -  Y-Visc  •  (3-14) 

Substituting  Eq.  (3.8)  and  (3.13)  into  Eg.  (3.7)  results  in 

°-£  '-    a^  • 

Choosing  a  freguency  of  100  Hz  for  the  worst  case  (greatest  viscous  and  thermal 
penetration  depths),  with  the  air  temperature  in  the  cavity  assumed  to  be  20°C, 
the  resulting  0surf  was  450.  Therefore,  viscous  and  thermal  wall  losses  are 
negligible. 
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Lastly,   the  loss  due  to  leakaqe  past  cavity  slideinq  surfaces  was 

considered.   Since  this  loss  could  not  be  directly  calculated,  a  measure  of 

its  significance  was  estimated  by  pressing  down  on  the  top  cover  of  the  cavity 

at  the  wall  edges  while  retaking  the   [100]  mode  resonance  curve  data  for  each 

sample.  The  empty  cavity  0  increased  from  147  to  252,  with  a  0.9  Hz  upshift  in 

the  center  frequency.   From  this  experiment,  a  value  of  Q  for  leakage  was 

estimated  to  be  approximately  280.   This  indicated  that  leakage  was  the  major 

source  of  loss  in  the  empty  cavity  by  at  least  a  factor  of  two.   The  leakage 

problem  could  not  be  corrected  without  making  major  changes  in  the  cavity 

design.  No  attempt  was  made  to  account  for  the  upshift  in  the  center  frequency 

that  occurred  while  pressinq  down  on  the  cavity  top  (the  frequency  upshift  was 

on  the  order  of  1  Hz  for  each  sample,  while  the  frequency  downshift  due  to 

introducing  the  fiberglass  was  on  the  average  of  10  to  12  Hz  drop  per  sample). 
4.  Determination  of  Rqc  from  fc  and  Q 

For  the  [100]  mode,  the  iterative  process  used  to  extract  the  DC  flow 
resistance  of  the  fiberglass  from  the  dispersion  relation  (Eq.  (2.51))  was 
started  by  choosing  Rqtj  equal  to  59,200  MKS  Rayl/m  (the  mean  flow  resistance 
determined  in  the  larqe  sample  holder  with  particle  velocity  normal  to  the 
air-fiberqlass  surface).  After  enterinq  the  fluid  and  cavity  parameters  (see 
Appendix  A  for  parameters  and  values  used),  the  roots  of  the  dispersion 
relation  for  ui  were  then  determined  for  each  sample  usinq  the  measured  values 
of  fc  and  Q  as  an  initial  starting  point.  The  resulting  roots  for  each  sample 
were  numerically  computed  until  the  sum  of  the  magnitudes  of  the  real  and 
imaginery  parts  of  the  dispersion  relation  (Eq.  2.51)  was  less  than  1  x  10~6. 
By  this  means,  calculated  values  of  Op  vs  fc  were  obtained  for  a  particular 
value  of  Roc* 

The  calculated  and  measured  values  of  0  vs  fc  were  compared  in  two 
ways.   Fiqure  3.16  shows  Oror  vs  fc*   Here  the  calculated  values  of  Qp  have 
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been  corrected  for  empty  cavity  losses  using  Eg.  (3.6).   Figure  3.17  shows  QF 
vs  fc.   Here  the  measured  values  of  Oj^jp  have  been  corrected  for  empty  cavity 

losses  using 

1      1       1 


0F    Otot    0e  • 

If  the  calculated  curve  plotted  above  the  measured  curve,  this 
indicated  too  high  a  value  of  R^  and  a  lower  value  of  R^  was  chosen.  If  the 
curve  plotted  below  the  measured  curve,  this  indicated  too  low  a  value  of  Rqtj 
and  a  larger  value  of  R^  was  chosen.  This  process  was  continued  until  the 
measured  curve  was  bounded  between  two  values  of  Rq£.  Once  bounded,  smaller 
increments  were  chosen  for  Rqtj  and  the  process  was  continued  until  a  best  fit 
(graphically)  was  obtained  between  the  calculated  and  measured  curves 
of  Of  or  Qt0T  vs  ^c*  For  tne  flOOJ  mode,  the  best  fit  value  for  R^  was 
15,000  P1KS  Rayl/m,  which  was  four  times  less  than  that  measured  experimentally. 

Two  possible  causes  were  investigated  to  explain  this  discrepency.  The 
first  possible  cause  is  the  10%  compression  of  the  samples  observed  in  the 
original  DC  flow  resistance  experiment  (large  sample  holder).  Assuming  the 
total  pressure  drop  across  the  test  sample  was  independent  of  the  compressed 
thickness,  the  extracted  value  of  DC  flow  resistance  per  unit  length  should  be 
corrected  from  59,200  MKS  Rayl/m  to  53,000  MKS  Rayl/m.  The  second  possible 
cause  is  that  the  flow  resistance  of  the  fiberglas  is  not  isotropic,  and  that 
it  is  greater  for  flow  perpendicular  to  the  sheet  faces.  This  might  be 
expected,  as  there  is  visual  evidence  that  the  fiberglass  sheets  have  been 
compressed.  In  the  first  flow  resistance  measurement,  the  fiberglass  sample 
was  oriented  so  that  the  fiberglass  sheet  surface  was  normal  to  the  flow 
velocity.  However,  for  the  [100]  mode,  the  particle  velocity  is  parallel  to 
the  fiberglass  sheet  surface.  Flow  resistance  tests  were  again  conducted, 
using  a  much  smaller  sample  holder,  with  samples  cut  so  that  the  anisotropy  in 
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RqC  could  be  measured.  As  indicated  earlier  in  this  chapter,  the  flow 
resistance  normal  to  the  sheet  surface  now  measured  approximately  three  times 
greater  than  that  parallel  to  the  sheet  surface.  The  value  measured  with  flow 
parallel  to  the  surface  was  15,200  MKS  Rayl/m.  With  this  value  of  &££,  the 
roots  of  the  dispersion  relation  were  recalculated,  Oitjt  and  ^F  calculated  and 
plotted.  These  values  fell  within  8%  of  measured  values  for  Qror  vs  ^c  anc^ 
within  12%  for  Qp  vs  fc  for  all  samples. 

The  method  of  extracting  Rqc  from  measured  values  of  fc  and  Q  was  then 
applied  to  the  [010]  mode,  which  has  particle  velocity  normal  to  the  sheet 
surface.  Starting  with  a  value  of  34,000  MKS  Rayl/m  for  R^rj  (the  value 
measured  in  the  small  sample  holder),  values  of  Q  vs  fc  were  calculated  for 
various  values  of  Rqq  and  were  compared  with  measured  values  of  Q  vs  fc.  The 
measured  values  of  Q  vs  fc  are  bounded  by  the  calculated  curves  for 
RDC  =  5'000  f1KS  Rayl/n  and  R^  =  20,000  MKS  Rayl/m,  and  are  consistent  with  a 
value  of  15,000  MKS  Rayl/m  for  R^f  the  sarne  as  £°r  the  [100]  mode.  This  is 
shown  in  Figures  3.18  and  3.19.  The  reason  for  the  discrepency  between  these 
bounding  values  and  the  measured  value  of  R^  for  flow  normal  to  the  fiberglass 
sheet  face  is  not  known.  It  may  be  due  to  mode  overlap  of  both  the  [110]  and 
[100]  modes  on  to  the  [010]  mode.  Figures  3.20  and  3.21  show  plots  of  the 
resonance  center  freguency  for  the  [100]  and  [010]  modes,  respectively,  as  a 
function  of  sample  thickness  for  the  same  values  of  Rq^  as  in  Figures  3.18 
and  3.19. 

Having  determined  that  the  best  fit  value  of  the  DC  flow  resistance  in 
the  porous  material  was  approximately  15,000  MKS  Rayl/m,  the  complex  wave 
number  and  characteristic  impedance  were  calcualted  using  Eq.  (2.56)  and  Eq. 
(2.57)  and  the  roots  of  the  dispersion  relation,  (Eq.  (2.51)).  Table  16  lists 
the  results  for  the  [100]  mode,  while  Table  17  lists  those  for  the  [010]  mode. 
The  data  for  both  modes  is  graphically  presented  in  Figures  3.22  through  3.26. 
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5.  Sensitivity  of  Roots  of  the  Dispersion  Relation  to  Variation  in 
Parameters 

In  order  to  determine  the  sensitivity  of  the  roots  of  the  dispersion 
relation  to  variations  in  fluid  and  porous  material  physical  properties,  the 
roots  were  again  calculated  for  the  [100]  mode  with  each  of  the  input 
parameters  changed  by  a  small  amount,  one  at  a  time,  and  the  characteristic 
impedance  and  complex  wave  number  were  recalculated. 

It  was  found  that  there  was  less  than  a  1%  change  in  k  and  Zc^ar 
when  porosity  was  decreased  to  0.95  (5%  change),  when  bulk  fluid  sound  speed 
was  changed  to  343  m/s  (0.6%  change  or  2  m/s),  and  when  bulk  fluid  density  was 
changed  to  1.23  Kg/m3  (1.6%  change  or  0.02  Kg/m3) .  However,  an  11%  change  in  k 
and  Zcnar  was  found  when  the  DC  flow  resistance  was  varied  by  ±  5,000  MKS 
Rayl/m  about  15,000  MKS  Rayl/m.  Because  the  uncertainty  in  the  measured  data 
is  so  much  greater,  the  sensitivity  of  the  best  fit  value  of  Rno  anci  hence  k 
and  Zcnar/  to  slight  uncertainties  in  other  physical  properties  is  negligible. 

6.  Modal  Pressure  Field  in  a  Rectangular  Cavity  with  Porous  Acoustic 
Material 

To  theoretically   determine  the  [100]   and   [010]   modal  pressure 

distribution  in  the  cavity  for  each  fiberglass  sample  reguires  the  values  of 

k^z  and  ^z     ^  determined  from  the  roots  of  the  dispersion  relation  for  to 

for  the  best  fit  value  of  Rq^.   Then,  using  Eq.  (2.63)  for  the  bulk  fluid  and 

Eq.  (2.64)  for  the  porous  material,  the  normalized  magnitude  and  phase  of  the 

pressure  at  each  grid  point  can  be  determined.   See  Appendix  B  for  a  more 

detailed  discussion.   For  the   [100]  mode,  Figure  3.27  shows  graphically  the 

results  for  sample  #1  (3/16")  and  Figure  3.28  shows  the  results  for  sample  #10 

(1  7/8").   For  the  [010]  mode,  Figure  3.31  shows  graphically  the  results  for 

sample  #1  (3/16")  and  Figure  3.32  shows  the  results  for  sample  #5  (11/16"). 
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Figure  3.22      |Zcharl/POc     vs     f     For  New  Method   ([100]    and    [010]    Modes) 
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Figure  3.23      4>zchar     vs     f     For  New  Method   ([100]    and    [010]    Modes) 
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Figure  3.24     k     vs     f     For  New  Method   ([100]    and    [010]    Modes) 
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Figure  3.25     a     vs     f     For  New  Method   ([100]    and    [010]    Modes) 
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Figure  3.26  cph  vs  f  For  New  Method  ([100]  and  [010]  Modes) 
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TABLE   16 

Zchar  AND  *  RESULTS   FDR    [100]    MODE 
(Rqq  =   15,200   NIKS  Rayl/m) 


Sample  # 

fc 
(Hz) 

k 
(Rad/m) 

a 

(Np/m) 

1  ^char 

|/Poc 

zchar 
(degrees) 

cph 
(m/s) 

1 

669.4  ±  0.1 

21.54 

14.53 

1.48 

-37.2 

195.2 

2 

665.4 

21.47 

14.47 

1.48 

-37.2 

194.8 

3 

657.7 

21.32 

14.35 

1.49 

-37.1 

193.8 

4 

649.8 

21.17 

14.21 

1.50 

-37.1 

192.9 

5 

637.6 

20.95 

13.98 

1.52 

-37.0 

191.2 

6 

625.3 

20.74 

13.71 

1.53 

-36.8 

189.4 

7 

617.1 

20.61 

13.52 

1.54 

-36.6 

188.1 

8 

601.5 

20.38 

13.10 

1.56 

-36.2 

185.4 

9 

587.4 

20.19 

12.66 

1.58 

-35.6 

182.8 

10 

569.5 

20.03 

11.96 

1.61 

-34.4 

178.6 

100 


TABLE    17 


Zchar  mD  %  RESULTS   FOR    [010]    MDDE 
(Rqq  =   15,000   MKS  Rayl/m) 


Sample  # 

fc 
(Hz) 

k 
(Rad/n) 

OL 

(Np/m) 

1 zchar 

/p0c 

^char 

(deqrees) 

cph 

(m/s) 

1 

865.0  ±  0.1 

25.76 

15.67 

1.30 

-34.8 

213.5 

2 

362.9 

25.72 

15.63 

1.30 

-34.8 

210.8 

3 

858.5 

25.64 

15.55 

1.31 

-34.8 

210.4 

4 

854.3 

25.57 

15.46 

1.31 

-34.8 

209.9 

5 

848.3 

25.47 

15.29 

1.32 

-34.7 

209.3 

101 


These  plots  present  the  node's  nodal  pressure  field  in  phasor  notation.  The 
naqnitude  is  represented  by  the  size  of  the  arrow  and  the  phase  angle  by  the 
anqular  rotation  of  the  arrow  in  degrees  from  the  horizontal  (i.e.,  zero 
degrees  of  phase  angle  is  a  phasor  pointing  horizontally  to  the  right) . 

Experimentally,  the  pressure  field  was  neasured  by  inserting  each 
sanple  into  the  cavity  and  exciting  the  desired  node  at  its  center  frequency 
for  that  sanple.  By  programing  a  HP-300  series  computer  to  take  data  from  the 
SR  530  Dual  Channel  Lock-In  Amplifier  (in  the  X,  Y  node),  the  in-phase  and 
guadrature  components  of  pressure  were  measured  at  each  grid  point.  The 
normalized  magnitude  and  phase  of  the  pressure  was  then  plotted.  For  the  [100] 
mode,  Figure  3.29  displays  the  results  for  sample  #1  (3/16")  and  Figure  3.30 
displays  the  results  for  sample  #10  (1  7/8").  The  results  for  the  [010]  mode 
are  displayed  in  Figures  3.33  and  3.34  for  sample  #1  (3/16")  and  sample  #5 
(11/16"),  respectively. 

In  comparing  the  theoretical  and  experimental  pressure  field  plots 
several  observations  can  be  made: 

(1)  The  effects  of  cavity  leakage  at  the  sliding  surfaces  is  evident 
when  the  experimental  plot  (Figure  3.29)  is  compared  to  the  theoretical 
plot  (Figure  3.27)  of  the  [100]  mode  for  sample  #1  (3/16")  by  noting 
that  the  pressure  phasors  alonq  the  boundaries  not  in  contact  with  the 
fiberqlass  strip  in  Figure  3.29  exibit  a  small  increase  in  phase  anqle 
(i.e.,  they  are  not  horizontal  as  they  would  be  if  no  losses  existed  at 
the  boundaries)  which  is  not  present  in  Fiqure  3.27.  This  effect  may 
also  be  seen  for  the  [010]  node  in  Fiqure  3.31  and  Fiqure  3.33. 

(2)  The  effects  of  the  porous  naterial  on  the  pressure  distribution  in 
the  bulk  fluid  can  be  seen  to  "qrow"  with  increasinq  sanple  thickness  by 
observinq  the  relative  increase  in  the  phase  anqle  away  from  the 
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air-f iberqlass  boundary  from  the  theoretical  plot  of   the  3/16"   (sanple 
ill)  fiberglass  strip  conpared  to  the  theoretical  nlot  of  the  1  7/8" 
(sample  #10)  fiberglass  strip.   This  can  be  seen  in  Figures  3.27  and 
3.28,  respectively.  This  effect  is  not  seen  in  the  theoretical  plots  of 
the  [010]  node  for  sample  #1  (Figure  3.31)   and  for  sanple   #5   (Figure 
3.32) 

(3)  The  effects  of  node  overlap  become  nore  evident  as  sample  thickness  is 
increased  by  comparing  the  measured  and  calculated  nodal  pressure  fields 
for  the  same  sample  thickness.  For  example,  compare  Figure  3.29  (sanple 
#1)  to  Figure  3.30  (sample  #10)  for  the  [100]  node  and  Figure  3.33 
(sanple  #1)  to  Figure  3.34  (sample  #5)  for  the  [010]  mode. 

(4)  The  pressure  distribution  at  the  air-fiberglass  boundary  shows  the 
effect  of  the  relative  size  and  sign  of  k\z  and  k2z  by  noting  that  for 
the  thickest  samples  the  pressure  phasor's  angle  increases  as  the 
interface  is  approached.  For  the  [100]  node  compare  the  theoretical  plot 
for  sample  #1  (Figure  3.27)  to  the  theoretical  plot  for  sample  #10. 
(Figure  3.28)  and  for  the  [010]  mode  compare  the  theoretical  plot  for 
sample  #1  (Figure  3.31)  to  the  theoretical  plot  for  sample  #5  (Figure 
3.32). 

The  primary  reason  for  comparing  the  calculated  and  measured  pressure 
distributions  was  to  verify  which  node  was  excited  and  to  check  its  purity.  The 
determination  of  whether  the  expected  node  had  been  excited  was  nade  by 
comparing  the  calculated  and  neasured  pressure  field  plots  for  the  thinnest 
fiberglass  strip  (3/16").  For  the  [100]  mode  the  nodal  plane  runs  normal  to  the 
length  direction  from  grid  point  (10,0)  to  grid  point  (10,15).  Travelling 
across  the  nodal  plane,  a  180  degree  phase  shift  in  the  pressure  phasor  is 
observed,  with  the  pressure  magnitude  reaching  a  maximum  at  the  cavity 
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boundaries  normal  to  the  length  direction.  Likewise,  for  the  [010]  mode  the 
nodal  plane  is  normal  to  the  width  direction  falling  half-way  between  the 
boundaries  parallel  to  the  length  direction.  The  pressure  phase  also  changes 
180  degrees  in  crossing  this  nodal  plane,  with  its  magnitude  reaching  a  maximum 
at  the  fiberglass  boundary  and  opposite  wall. 

The  purity  of  these  modes  is  retained  only  with  the  thinner  samples  and 
becomes  contaminated  as  thicker  samples  are  introduced.  This  is  clearly  seen  by 
comparing  Figures  3.29  and  3.30. 
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IV.   COMPARISON  OF  RESULTS  AND  CONCLUSION 

A.   COMPARISON  OF  THE  COMPLEX  W/E  NUMBER 

The  complex  wave  number  may  be  written  in  terms  of  its  real  and  imaqinery 
parts  as  the  wave  number  k  and  the  absorption  coefficient  a  ,  respectively. 
Figure  4.1  shows  the  wave  number  k  as  a  function  of  frequency  for  each  method. 
Values  of  k  determined  by  the  new  method  fall  between  those  found  by  Yaniv's 
method  (the  upper  bound)  and  by  Beranek's  Material  Properties  method  in  the 
isothermal  limit  (the  lower  bound).  For  example,  at  600  Hz  the  value  of  the  wave 
number  determined  by  the  new  method  was  within  6%  of  the  value  determined  by 
Yaniv's  method.  The  values  of  k  from  the  probe  tube  experiment  are  about  twice 
those  determined  by  the  new  method.  This  is  to  be  expected,  since  the  particle 
velocity  in  the  probe  tube  experiment  was  normal  to  the  air-fiberglass  boundary, 
i.e.,  was  in  the  direction  of  hiqher  measured  flow  resistance.  (Note  that 
although  the  [010]  mode  also  has  particle  velocity  normal  to  the  fiberglass 
surface,  the  best  fit  value  of  Rqq  (15,000  MKS  Rayl/m)  was  used  to  calculate  both 
the  complex  wave  number  and  characteristic  impedance). 

The  abosorption  coefficient  as  a  function  of  frequency  is  displayed  in  Fiqure 
4.2  for  each  method.  The  values  determined  by  the  new  method  fall  about  those 
obtained  by  Beranek's  Material  Properties  method  (in  the  isothermal  limit).  For 
example,  at  600  Hz  the  value  of  the  absorption  coefficient  determined  from  the 
new  method  falls  with  17%  of  the  value  determined  by  Yaniv's  method.  Again,  the 
values  determined  usinq  the  probe  tube  method  were  about  twice  those  measured  by 
the  new  method.  The  cause  of  the  dip  in  the  values  of  the  absorption  coefficient 
determined  by  Yaniv's  method  is  not  known. 

Fiqure  4.3  displays  the  phase  speed  in  the  porous  material  as  a  function  of 
frequency  for  each  method.   Values  determined  by  the  new  method  are  bounded  by 
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those  found  by  Yaniv's  method  and  by  Beranek's  Material  Properties  method  (in 
the  isothermal  limit).  For  example,  at  600  Hz  the  phase  speed  determined  by  the 
new  method  is  within  7%  of  the  value  determined  by  Yaniv's  method.  The  values 
of  phase  speed  determined  by  the  probe  tube  method  are  approximately  one-half  of 
those  found  by  the  new  method  for  the  same  reason  discussed  above. 

B.  COMPARISON  OF  THE  CHARACTERISTIC  IMPEDANCE 

The  normalized  characteristic  impedance  as  a  function  of  frequency  for  all 
methods  is  displayed  in  Figure  4.4.  Values  determined  by  the  new  method  are 
bounded  by  those  from  the  adiabatic  and  the  isothermal  material  properties 
method.  Values  from  the  new  method  plot  almost  directly  on  top  of  the  Biot 
theory  in  the  diffusion  limit.  For  example,  at  600  Hz  the  value  determined  by 
the  new  method  fell  within  26%  of  that  determined  by  Yaniv's  method  and  within 
20%  of  the  impedance  tube  value. 

Figure  4.5  displays  the  impedance  phase  as  a  function  of  frequency  for  each 
method.  The  values  determined  by  the  new  method  fall  about  the  values 
determined  by  the  material  properties  method.  The  low  frequency  values  (<  200 
Hz)  determined  by  the  impedance  tube  method  are  thought  to  be  in  error  due  to 
the  difficulties  in  measuring  the  distance  to  the  first  node  in  that  method.  It 
is  not  known  why  the  values  of  the  phase  extracted  using  the  new  method  from  the 
[100]  mode  data  increase  for  decreasing  freguency,  since  in  the  low  frequency 
limit,  the  phase  should  approach  -45  degrees  (pure  diffusion  of  sound  pressure). 

C.  CONCLUSION 

The  ability  to  determine  a  material's  DC  flow  resistance  from  parameters  of 

the  resonance  curve  in  a  plane  wave  resonator  was  verified.   Additionally,  the 

values  of  both  the  complex  wave  number  and  characteristic  impedance  determined 

from  the  new  method  compared  well  with  classical  methods;  the  values  obtained  by 

the  various  methods  agree  with  each  other  to  within  26%. 
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APPENDIX  A 

METHOD  OF  SOLUTION  OF  THE  DISPERSION  RELATION  \ND  ITS  IMPLEMENTATION 

The  mehtod  used  to  find  the  roots  of  the  dispersion  relation  (Eq.  (2.51)) 
was  Ward's  nethod  [Ref .  19,  pp.  169-170] .  This  method  is  not  the  only  method 
available  to  find  the  roots  of  a  complex  transcendental  equation  but  is  easy  to 
understand  and  to  implement,  as  no  deriviatives  are  needed. 

1.   WARD'S  METHOD 

If  we  let  the  polynomial  f(z)  be  a  function  of  the  complex  variable 
z  =  x  +  iy  such  that 

f(z)  =  u(x,y)  +  iv(x,y)  ,  (A. 1.1) 

it  can  be  shown  that  f(z)  has  a  zero  a  =  b  +  id  if  and  only  if  u(b,d)  =  v(b,d)  = 
0.  This  pair  of  simultaneous  equations  in  u(x,y)  and  v(x,y)  may  be 
conveniently  solved  in  two  ways.  First,  one  can  seek  to  minimize 

u2(x,y)  +  v2(x,y)  (A. 1.2) 

or,  secondly,  on  can  minimize 

w(x,y)  =  |u(x,y)|  +  |v(x,y)|.  (A. 1.3) 

Since  no  derivitives  are  employed,  this  method  seems  well  suited  for  findinq 
multiple  roots. 

Ward's  technique  is  based  on  the  knowledqe  that  in  every  neiqhborhood  of  a 
point  (x,y)  such  that  w(x,y)  *  0,  there  exists  a  point  w(x',y')  such  that  w(x', 
y')  <  w(x,y).  At  a  qiven  step  in  the  process,  which  is  iterative,  the  current 
value  of  w(x,y)  is  compared  with  the  vlues  w(x+h,y),  w(x-h,y),  w(x,y+h)  and 
w(x,y-h)  until  a  smaller  value  is  found  (if  possible).  The  coordinates  yieldinq 
this  smaller  value  of  x  and  y  are  used  in  the  next  iterative  step.  Should  these 
four  points  yield  no  decrease  in  w  =  |u|  +  |v|,  the  increment  h  is  decreased 
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by  some  amount  and  the  process  is  repeated.  As  with  other  methods  which  can  sho\ 
convergence  within  certain  limits,  there  is  no  guarantee  of  convergence 
with  this  method. 

2.   IMPLEMENTING  WARD' S  METHOD 

In  order  to  apply  Ward's  method  to  find  the  roots  of  the  dispersion 
relation,  it  must  be  separated  into  its  real  and  imaginary  parts.  To  do  so  note 
that  the  parameters  k]_z  and  k2z,  and  a(u)  are  functions  of  w  ,  they  may  be 
separated  into  their  respective  real  and  imaginary  components  (these  three  terms 
also  have  complex  parameters  buried  within  them  which  are  freguency  dependent). 
Although  tedious,  it  is  necessary  with  this  method. 

Once  each  parameter  has  been  reduced  to  its  simplest  complex  form  (real 
and  imaginary  portions),  these  can  be  inserted  back  into  the  dispersion  eguation 
and  by  use  of  complex  arithmetic,  the  dispersion  relation  may  be  recast  in  the 
form 

Re  {Dispersion  Relation}  +  jIm{Disperson  Relation}  =  0 
and  Ward's  Method  of  Eg.  (A. 1.3)  may  be  now  implemented  in  a  computer  algorithm. 

The  numerical  analysis  program  written  using  Ward's  Method  reguired  the 
following  inputs  for  each  sample  thickness: 

( 1)  Measured  Quantities 

a)  fc  (measured  in  the  cavity) 

b)  0   (measured  in  the  cavity) 

(2)  Cavity  Parameters 

a)  Length  of  the  cavity,  L.  (in  x-direction) 

b)  Width  of  porous  material,  d2.   (sample  thickness) 

c)  Width  of  remaining  bulk  fluid,  d]_. 

d)  Mode  number  in  the  long  direction.   {n  =  0  for  the   [010]  mode, 
n  =  1  for  the  [100]  mode} 
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(3)  Fluid  and  Porous  Solid  Parameters 

a)  Adiabatic  fluid  sound  speed,  Cf.   (<§  23°C,  345.2  m/s) 

b)  Fluid  bulk  density,  pf.   ((§  20°C,  1.21  Kq/m3) 

c)  Adiabatic  ratio  of  specific  heat  of  fluid,  y.  (<a  20°C,  1.402) 

d)  Prandtl  number,  Pr.   (@  20°C,  0.733) 

e)  Porosity,  P.  (  *  1.0) 

(4)  Program  Parameters 

a)  Angular  frequency  increment,  Aw. 

b)  Damping  freguency  increment,  AS. 

c)  An  ending  criterion  for  w(  oj,  g)  . 

(5)  Initial  Estimate  of  Rpr 

a)  If  no  a  priori  knowledge  of  material's  flow  resistance,   then 
assume  very  hiogh  or  low  value  as  the  starting  point. 

b)  Perform  a  DC  flow  resistance  experiment  to  extract  a  mean  value 
for  Rqc- 

Once  these  quantities  are  entered,  the  program  calculates  its  first  w( u, g)  based 
on  the  initial  inputs  of  fc  and  Q.  This  value  of  w(w,  g)  will  be  compared  to  the 
four  incremental  values  of  w(  co  ±  Ato,  g  ±  AS),  that  is,  it  compares  the 
initial  w(  oo,  S)  with 

w(  orl-Aoj,  8) 

w(  w-Au),  S) 

w(w,  8+AB) 

w(w,  S-AS) 
Which  ever  is  the  lowest,  those  coordinates  in  w, S  space  are  chosen  as  the  new 
values  and  the  process  is  continued  until  w(a),  g)   is  less  than  the  ending 
criterion.  The  output  converts  the  roots  of  u>  and  g  to  fc  and  Q. 
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This  process  is  continued  for  each  sample  thickness,  and  the  resulting 
values  can  be  plotted  either  as  0  vs  fc  or  as  fc  vs  sample  thickness  and 
compared  with  the  respective  curve  after  loss  corrections  to  Q  are  made.  An 
evaluation  is  made  as  to  whether  Rqq  *s  to°  high  or  too  low,  the  value  of  R^-. 
changed  and  the  process  repeated  for  this  new  value  of  Rqc  for  all  samples. 
This  iterative  process  is  completed  when  the  roots  are  found  for  a  value  of  Rqq 
which  provides  the  best  fit  curve  through  the  measured  data.  With  this  value  of 
R££  and  the  roots  for  each  sample  thickness  determined,  the  characteristic 
impedance  and  complex  wave  number  may  be  calculated  by  use  of  Eq.  (2.56)  and  Eq. 
(2.57). 
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APPENDIX  B 

DETERMI NATION  OF  TOE  THEORETICAL  MODAL  PRESSURE  DISTRIBUTION  IN  A  PLANE 
WAVE  RESONATOR  WITH  AN  ACOUSTICALLY  ABSORBING  MATERIAL 


The  equation  that  governs  the  pressure  distributuion  in  the  bulk  fluid 
(region  1)  is: 


PI 


-joopfcos[kxx]  Jig  cos[klz(z-d1)] 


jwt 


(B.l.l) 


The  equation  that  governs  the  pressure  distribution  in  the  porous  material 
(region  2)  is,  after  substituting  Eq.  (2.62)  for  4>20: 


'P2 


cos[k.    d.  ]  ^ 

■j^pfcos[kxx]  ?10  cosRfVl"  cos^2z(z+d2)] 

2z  2 


jut 


(B.1.2) 


With  the  value  of  Rq^  known,  along  with  the  roots  of  the  dispersion  relation  for 
a  given  sample  thickness,  the  pressure  distribution  in  region  1  and  region  2 
can  be  calculated. 

In  the  experimental  measurement  of  the  modal  pressure  field,  the  quadrature 
component  of  the  pressure  was  zeroed  at  the  reference  position  of  (0,0,1)  to  set 
a  reference  point  such  that  both  theory  and  experiment  could  be  compared.  The 
0,0  indicates  the  grid  point  (0,0)  and  the  1  indicates  that  the  receiver  probe 
tube  was  located  in  the  ceiling  of  the  cavity.  To  set  the  reference  at  x  =  0,  z 
=  di  such  that  the  phase  of  the  pressure  distribution  at  this  point  is  zero, 
requires  that 

4>10  =  1  /-ATN(-u)/8)   .  (B.1.3) 

To  calculate  the  theoretical  modal  pressure  distribution,  start  at  the 
reference  location  of  x  =0  and  z  =  d^  and  incrementally  move  away  from  this 
position  in  increments  of  Ax  and  Az  such  that  the  theoretical  grid  points  match 
those  existing  in  the  actual  cavity. 
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